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ABSTRACT: Random matrix theory is a powerful way to describe universal correlations of 
eigenvalues of complex systems. It also may serve as a schematic model for disorder in quantum 
systems. In this review, we discuss both types of applications of chiral random matrix theory 
to the QCD partition function. We show that constraints imposed by chiral symmetry and 
its spontaneous breaking determine the structure of low-energy effective partition functions for 
the Dirac spectrum. We thus derive exact results for the low-lying eigenvalues of the QCD 
Dirac operator. We argue that the statistical properties of these eigenvalues are universal and 
can be described by a random matrix theory with the global symmetries of the QCD partition 
function. The total number of such eigenvalues increases with the square root of the Euclidean 
four-volume. The spectral density for larger eigenvalues (but still well below a typical hadronic 
mass scale) also follows from the same low-energy effective partition function. The validity of 
the random matrix approach has been confirmed by many lattice QCD simulations in a wide 
parameter range. Stimulated by the success of the chiral random matrix theory in the description 
of universal properties of the Dirac eigenvalues, the random matrix model is extended to nonzero 
temperature and chemical potential. In this way we obtain qualitative results for the QCD phase 
diagram and the spectrum of the QCD Dirac operator. We discuss the nature of the quenched 
approximation and analyze quenched Dirac spectra at nonzero baryon density in terms of an 
effective partition function. Relations with other fields are also discussed. 
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1 INTRODUCTION 

Well before the advent of quantum chromodynamics (QCD), the theory of the 
strong force, it was realized that the essential ingredients of the hadronic world at 
low energies are chiral symmetry and its spontaneous breaking, (see e.g. Refs. [jl], 
H). Mainly through lattice QCD simulations, it has become well established 
that chiral symmetry breaking by the vacuum state of QCD is a nonperturbative 
phenomenon that results from the interaction of many microscopic degrees of 
freedom. We argue in this review that the complexity of the QCD vacuum leads 
to a low-energy description that is completely dictated by the global symmetries 
of QCD. This interpretation of Goldstone's theorem provides a natural duality 
between a strongly interacting fundamental theory and a weakly interacting low- 
energy effective theory. 

1.1 QCD and Chiral Symmetry 

We illustrate the concept of spontaneous symmetry breaking using the simpler 
example of a classical spin system with two rotational degrees of freedom. The 
Hamiltonian of this system has a certain symmetry: It is invariant under rota- 
tions. In mathematical language, the symmetry group is G = 0(3). However, 
at low temperatures, the ground state of the system does not exhibit this sym- 
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metry. In a small external magnetic field, which breaks the rotational invariance 
explicitly, the spins will polarize in the direction of the magnetic field. In the 
thermodynamic limit, the spins remain polarized even if the magnetic field is 
switched off completely. This phenomenon is known as spontaneous magnetiza- 
tion. The ground state is no longer invariant under 0(3) rotations but only under 
0(2) rotations in the plane perpendicular to the spontaneous magnetization. In 
mathematical language, the full symmetry group G = 0(3) is spontaneously bro- 
ken to a smaller symmetry group H = 0(2). The spontaneously broken phase 
is characterized by low-energy excitations in the form of spin waves in the plane 
perpendicular to the spontaneous magnetization. This is a consequence of a gen- 
eral theorem known as Goldstone's theorem ||, which tells us that spontaneous 
breaking of a continuous symmetry leads to low- lying excitations, the Goldstone 
modes, with a mass that vanishes in the absence of a symmetry-breaking field. 

The Goldstone modes are given by the fluctuations in the plane perpendicular 
to the direction of the spontaneous magnetization. Thus, the spin system has 
two Goldstone modes. In general, spontaneous symmetry breaking in a system 
of spins with n components is associated with n — 1 Goldstone modes. This 
number is also equal to the number of generators of the coset G/H [the number 
of generators of 0(n) is n(n — l)/2]. 

A spontaneously broken symmetry is characterized by an order parameter, 
which in this case is the spontaneous magnetization. At nonzero temperature, 
the alignment of the spins is counteracted by their thermal motion, and above 
a critical temperature (the Curie temperature) the spontaneous magnetization 
vanishes. 

Let us now consider the hadronic world and interpret the particle spectrum 
in terms of the concepts discussed above. We will look for the simplest theory 
consistent with the following two empirical facts: (a) there are three bosonic 
particles, the pions, that are much lighter than all other hadrons and (b) the 
proton and the neutron have almost the same mass. Fact (a) implies that there 
are three Goldstone bosons associated with a spontaneously broken symmetry. 
Assume that the underlying theory is based on an n-component field without 
preferred directions, i.e. the theory is invariant under 0(n) transformations of 
the fields. Spontaneous symmetry breaking means that the ground state has 
a preferred direction, leaving n — 1 directions for the Goldstone modes. Since 
there are three pions, this suggests that n = 4 and that the ground state of 
the theory should be symmetric under 0(3). The 0(3) symmetry is the familiar 
isospin symmetry, which results in the near equality of the masses of the pions 
and of the mass of the proton and the neutron, fact (6). As was first conjectured 
by Gell-Mann and Levy ^j, the theory of the strong force is based on an 0(4) 
invariance spontaneously broken to 0(3) with nucleons transforming according 
to SU(2). This is the familiar linear a- model. Since the pions are not completely 
massless, a small 0(4) symmetry-breaking mass term should also be present in 
the theory. 

The o"-model is a phenomenological model for the interactions of pions and nu- 
cleons. It has been very successful in explaining many previously known empirical 
relations. Our aim, however, is to understand chiral symmetry in terms of QCD, 
the fundamental theory of the strong interactions. Quantum chromodynamics 
is a gauge theory of quarks that come in six different flavors and three different 
colors. The gauge field interaction is according to the non-Abelian SU(3) color 
group. Two of the six quarks are nearly massless, and at low energies QCD is 
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well approximated by a theory with only the two lightest quarks. They are mixed 
by the SU(2) isospin symmetry group. This symmetry is exact for degenerate 
quark masses. For massless quarks, there is an additional symmetry. The helicity 
of a particle is a good quantum number, and the right-handed and left-handed 
quarks can be rotated independently. The isospin or vector symmetry rotates 
both chiralities in the same way whereas the axial SU(2) group rotates them 
in the opposite direction. This explains that the chiral symmetry group in the 
massless case is G = SU(2) x SU(2), which is isomorphic to 0(4). The mass 
term breaks this symmetry explicitly and thus plays the role of the magnetic 
field in the spin system. Even in the massless case, however, the vacuum state 
of QCD is characterized by a nonzero expectation value of the chiral condensate, 
which, like the mass term, mixes right-handed and left-handed quarks. Thus, in 
the vacuum, the axial part of the symmetry group G is broken spontaneously. 
The ground state is not unique, and the degenerate states are connected by the 
broken group G. The degeneracy can be lifted by means of a small mass term, 
and in the thermodynamic limit the system will be frozen in this direction, i.e. 
the direction of the QCD vacuum state is determined by the mass term. This is 
exactly the same situation as in the spin system discussed above. The Goldstone 
modes analogous to spin waves are the pions, with a mass that is well below the 
typical hadronic mass scale of about 1 GeV. As in the spin system, we expect that 
the expectation value of the chiral condensate will become zero above a critical 
temperature. This phenomenon is known as the restoration of chiral symmetry. 
The chirally symmetric phase probably existed in the early universe and may be 
produced in relativistic heavy-ion collisions. 

Although QCD is the consistent theory of the strong interactions, many ques- 
tions remain unanswered. For example, why have its constituents never been 
observed in nature? This phenomenon is known as confinement and means that 
all physical states are color singlets. QCD is best understood at high energies, 
where, because of asymptotic freedom, quarks and gluons become weakly inter- 
acting and perturbative calculations are possible. At low energies, on the other 
hand, it is necessary to rely on nonperturbative approaches. One approach is to 
study the QCD partition function by means of Monte Carlo simulations of a dis- 
cretized version of the QCD action. This approach has been very fruitful, and a 
great deal of our understanding of low-energy QCD is based on such calculations 
(for a recent review, see Ref. ||]). The drawback is that large-scale simulations 
do not necessarily provide a simple picture of the relevant degrees of freedom. 
Therefore, it is often advantageous to study QCD by means of effective models 
or theories. One example is the instanton liquid vacuum ||, which is a model for 
an ensemble of relevant gauge field configurations. Effective low-energy theories 
are a second example. We already mentioned the cr-model for pions and nuclc- 
ons. However, because of the spontaneous breaking of chiral symmetry and the 
existence of a mass gap in QCD, one can do better. Based on chiral symmetry 
one can formulate an exact low-energy theory for the Goldstone modes. This 
nonlinear cr-model is the basis for a systematic low-energy expansion || [?] that 
is discussed in detail in Sees. L3 and |I[ 



The QCD partition function can be written as a Euclidean path integral that 
can be expressed as the expectation value of the fermion determinant, 

Z = (j[det(V + m f )) . (1) 
/ 
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Here the average is over all gauge fields weighted by the Euclidean Yang-Mills 
action, the product is over quark flavors of mass m/, and T> is the Dirac op- 
erator, which we introduce in great detail in Sec. ^. The fermion determinant 
can be expressed as a product over the eigenvalues i\ n of the Dirac operator. 
Therefore, we may also interpret the average as an average over the eigenvalues 
with probability distribution determined by the gauge field dynamics. We have 
argued that spontaneous breaking of chiral symmetry means that a small quark 
mass leads to a macroscopic realignment of the QCD vacuum. It is clear from the 
QCD partition function that this is only possible if there is an accumulation of 
Dirac eigenvalues near zero. Otherwise, a small mass term would be completely 
dominated by the much larger eigenvalues in the factors of (i\ n + nif). For a free 
Dirac operator in a four-dimensional box, the eigenvalue density is proportional 
to A 3 and thus vanishes near zero. Therefore, the small eigenvalues must be due 
to interactions mediated by the gauge fields. 

There are two possibilities. First, the eigenvalues may originate from the ex- 
actly zero Dirac eigenvalues in the field of an instanton. At a nonzero density of 
the liquid of instantons and anti-instantons, the zero eigenvalues are distributed 
over a band because of interactions that lift the degeneracy of the eigenvalues 
|||. The eigenvalue repulsion results in a nonzero eigenvalue density near zero. 
Second, the eigenvalues may originate from the bulk of the Dirac spectrum. As 
is the case for any interaction in quantum mechanics, the interactions mediated 
by the gauge fields lead to a repulsion of the eigenvalues. For increasing inter- 
action strength, it is then advantageous for the eigenvalues to move to a region 
with a low eigenvalue density. Both mechanisms of spontaneous chiral symmetry 
breaking rely on the repulsion of eigenvalues in interacting systems. This topic 
has been investigated in detail by means of random matrix theory, discussed in 
the next subsection. 

As mentioned above, we expect chiral symmetry to be restored at high temper- 
atures Q. This is confirmed by lattice QCD simulations, which show a critical 
temperature of about 150 MeV ||10|| . However, there is another direction in which 
the QCD partition function can be explored, namely at nonzero baryon number 
density. In this domain, rigorous results are available only at infinite baryon 
number density |yj. This raises the question to what extent these results can 
be extrapolated to physically interesting densities. Because of the complex phase 
of the fermion determinant, Monte Carlo simulations are not possible in this 
case. Therefore, it is not surprising that many of the recent developments in 
QCD at nonzero density are based on the analysis of effective models ]l2|, [D| . 
The picture that has emerged for two massless flavors is that a first-order chiral 
phase transition occurs at zero temperature along the chemical potential axis. 
Renormalization-group arguments and lattice QCD simulations indicate that for 
two massless flavors, a second-order phase transition occurs at zero chemical po- 
tential along the temperature axis. The expectation is that a line of first-order 
transitions and a line of second-order transitions will extend in the \xT plane and 
will join at the tricritical point, as indicated in Fig. 111. A color-superconducting 
phase is conjectured to exist at higher densities [12, [13, 11] but the discussion of 



this phase is beyond the scope of this review. One of the questions we address 
is the robustness of this picture based on the dynamics of the eigenvalues of the 
QCD Dirac operator. 
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Figure 1: A minimal phase diagram for QCD with two massless flavors at nonzero 
temperature and density. 



1.2 Random Matrix Theory 

Random matrix theory (RMT) first appeared in the mathematical literature in 
1928 [jl4| and was first applied to physics in the context of nuclear resonances by 
Wigner almost 50 years ago |15[| . At that time, theoretical approaches such as 
the shell model had been very successful in describing the low-lying excitations 
of complex nuclei. However, highly excited resonances, which can be observed 
experimentally by neutron scattering, could not be described by the microscopic 
theory. The problem is generic: For any complex quantum system containing 
many degrees of freedom with complicated dynamics, it is very hard, if not im- 
possible, to obtain exact results for the energy levels far above the ground state 
of the system. 

Having acknowledged that the highly excited states cannot be predicted indi- 
vidually, one can ask whether the experimental data have some generic statistical 
features that can be described theoretically. This is where RMT comes in. Every 
quantum system is described by a Hamilton operator that can be expressed in 
matrix form. For a complex quantum system such as a large nucleus, this matrix 
is very complicated, and we may not even know its details. In this case, one ap- 
proach is to assume that all interactions that are consistent with the symmetries 
of the system are equally likely. This means replacing the elements of the Hamil- 
ton matrix by random numbers that are uncorrelated and distributed according 
to the same probability distribution. To obtain definite results, observables such 
as the level density must then be averaged over the random matrix elements. 
This defines a statistical theory of energy levels, which is mathematically known 
as RMT. An important point is that the random matrix must have the same sym- 
metries as the original Hamilton matrix. A collection of early papers on RMT 
can be found in the book by Porter [16], and the standard reference on RMT is 
the book by Mehta 0. 

Can such an enormous simplification of the real problem describe empirical 
data? A similar puzzle occurred in the theory of critical phenomena, in which 
the critical behavior does not depend on the detailed dynamics of the theory. The 
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reason for this simplification is the appearance of a length scale, the correlation 
length, that diverges at the critical point. Because of the corresponding sepa- 
ration of scales, it is possible to integrate out the short-wavelength fluctuations 
and renormalize the theory to a fixed-point theory that does not depend on the 
details of the initial theory. What is the separation of scales that takes place in 
quantum spectra? The two basic scales are the average level spacing and the scale 
of the variation of the average level spacing. The equivalent of the correlation 
length is the inverse average level spacing, which diverges in the thermodynamic 
limit or the semiclassical limit. We thus expect that spectral properties on the 
scale of the average level spacing do not depend on the details of the underlying 
dynamics. They are universal. 

Universal properties can be studied in the simplest model of a given universality 
class and, in this way, exact analytical results for the correlation functions can 
be obtained. In the case of spectral correlations the simplest models of the 
universality classes are the Gaussian RMTs, and in the past three decades many 
exact results have been derived for these models. Unfortunately, in most cases 
it is not possible to prove whether a given theory belongs to one of the RMT 
universality classes. Therefore, random matrix predictions are usually verified by 
comparisons with empirical data. 

Since Wigner's original proposal, universal quantities have been identified and 
computed in a variety of fields including nuclear physics, atomic and molecular 
physics, disordered mesoscopic systems, quantum systems with classically chaotic 
analogs, two-dimensional quantum gravity, conformal field theory, and QCD. A 
recent comprehensive review of the applications of RMT can be found in Ref. [^] . 
RMT is now an independent subfield of mathematical physics. It provides a 
unifying description of universal statistical features of many different quantum 
systems and is applicable whenever a system is sufficiently complex. 

Let us raise an interesting point here. The eigenvalues of the Hamilton matrix 
for a given quantum system are the quantities of primary interest. Instead of 
going through RMT to compute the eigenvalues, an alternative — and much 
simpler — approach might be to postulate random eigenvalues. It turns out 
that this does not describe empirical data, at least not if the random eigenvalues 
are uncorrelated. Therefore, in addition, one would have to postulate how the 
eigenvalues are correlated, and it is not at all clear how to do this. In RMT, 
on the other hand, one starts with uncorrelated random matrix elements. At 
the end of the calculation, one finds that the resulting eigenvalues are strongly 
correlated in precisely the right way to describe the data. 

In most applications, RMT is used to describe the statistics of energy levels, 
i.e. of the eigenvalues of the Hamilton operator. For the Euclidean QCD partition 
function, it is more natural to construct a random matrix model for the Dirac 
operator. As shown in the previous subsection, the spectrum of this operator is 
intimately related to the phenomenon of chiral symmetry breaking. This estab- 
lishes the connection between RMT and chiral symmetry announced in the title 
of this review. We use random matrix methods to study the Dirac spectrum and 
its implications for chiral symmetry breaking. 

It is important to note that RMT can only provide an exact description of uni- 
versal quantities. It cannot be used for the calculation of nonuniversal observables 
such as the average level spacing. Such behavior is well known in statistical me- 
chanics where, for example, the Ising model describes the critical exponents of the 
liquid-gas phase transition but does not give the critical temperature. In the case 
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of the distribution of the eigenvalues of a system, the global spectral density is 
not described by RMT. For example, for the Gaussian random matrix ensembles 
the average spectral density is a semicircle, whereas in real systems it typically 
increases strongly with excitation energy. In contrast, universal quantities do not 
depend on the details of the dynamics of the system. In RMT, they do not depend 
on the distribution of the random matrix elements. It is crucial to distinguish the 
universal quantities from the model-dependent ones. In the first five sections of 
this review, we mainly concentrate on the universal, model-independent features 
that yield exact quantitative results. However, it is sometimes useful to construct 
random matrix models to obtain a qualitative description of the physics in a dis- 
ordered system. As we show in Sec. ^, schematic random matrix models for QCD 
at nonzero temperature and density can yield important qualitative insights into 
problems that are otherwise difficult to tackle. 



1.3 Effective Low-Energy Theories and Chiral Random Matrix Theory 

Let us now make the connection between QCD, effective low-energy theories for 
QCD, and chiral RMT. We concentrate here on the big picture, postponing a 
more detailed discussion to Sees. ||] and f|. 

At low energies, quarks and gluons are confined in hadrons, i.e. the particles 
we observe in nature are composite objects. Instead of attempting to describe 
the results of low-energy scattering experiments in terms of quarks and gluons, 
it is often simpler to use an effective theory whose elementary degrees of freedom 
are the lightest particles of the theory. As mentioned earlier, in QCD the low- 
lying degrees of freedom are the pions resulting from the spontaneous breaking 
of chiral symmetry. Therefore, an essential ingredient in the construction of the 
effective low-energy theory is the requirement that it correctly incorporates the 
chiral symmetries of the original theory. Since the up- and down-quark are not 
completely massless in nature, the pions are not massless but have a small mass of 
about 140 MeV. They are much lighter than the lightest non-Goldstone particles, 
such as the p-meson or the nucleons, which have a mass of about A ~ 1 GeV. 
This means that at sufficiently low energies, the QCD partition function is well 
approximated by the partition function of an effective low-energy theory involving 
only pions. 

We now consider QCD in a finite Euclidean volume V4 = L 4 . The partition 
function is then dominated by the pions if 

j«L. (2) 

This statement follows from Eq. (|j) below by comparing the contribution of the 
pion, exp(— m n L), to that of a heavier particle, exp(— AL). The interactions 
of the pions are described by an effective chiral Lagrangian that will be given 
in Sec. £§ The fields in this La grangian can be separated into zero-momentum 
modes (constant fields) and nonzero-momentum modes. It was realized by Gasser 
and Leutwyler |19[ | that there exists a kinematic regime where the fluctuations of 
the zero-momentum modes dominate the fluctuations of the nonzero-momentum 
modes. This regime is given by the condition 

L < — , (3) 
rriir 
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where m n is the pion mass. Intuitively, this means that the wavelength of the 
pion is much larger than the linear extent of the box. Thus, the pion field does 
not vary appreciably over the size of the box, so the derivative terms are small. 
(These statements are quantified in Sec. 0.) We need consider only the zero- 
momentum modes in the regime of Eq. (|!|). For constant fields, the spacetime 
integral in the action can be replaced by the four-volume, and, effectively, we 
only have to deal with a much simpler zero-dimensional theory. However, the 
global symmetries remain important. 

We are now ready to make the connection to chiral RMT. A random matrix, as 
introduced in the previous subsection, contains independently distributed random 
variables. Each element of the matrix has the same average size. On the other 
hand, the matrix elements of the Dirac operator contain specific correlations due 
to the gauge fields and the kinetic terms, which are not included in the RMT. 
Only the global symmetries of QCD are included in the RMT. Since the random 
matrix elements do not have any spacetime dependence, we expect that in the 
domain of Eq. (|3|), which is dominated by constant fields, the RMT reproduces 
the mass dependence of the finite- volume partition function. 

To summarize, the three ingredients in the construction of the finite-volume 
effective theory from QCD in the regime of Eq. @ are the following: (a) the 
global symmetries of QCD, (6) the spontaneous breaking of chiral symmetry, and 
(c) the fact that the partition function is dominated by the constant Goldstone 
fields. There is an almost one-to-one correspondence to the properties of the chiral 
random matrix model (2(J : (a) the Dirac matrix is replaced by a random matrix 
with the same global symmetries, (b) chiral symmetry is broken spontaneously in 
the limit of infinitely large random matrices, and (c) the random matrix elements 
do not have any spacetime dependence. 

2 QCD PARTITION FUNCTION AND DIRAC SPECTRUM 

After the introductory remarks on our basic philosophy in Sec. |], we now present 
a detailed discussion of the properties of the QCD partition function and the 
Dirac operator. 

2.1 Basic Definitions 

The QCD partition function is defined as 

Z QCD = Tre -/3H ; (4) 

where (3 is the inverse temperature and H is the QCD Hamiltonian in a box 
of volume V3 = L 3 . It can be rewritten as a (suitably regularized) functional 
integral in Euclidean space, 

r N t 

Z QCB = / DA ^ -q det( p + m ^ e -s YM ; ( 5 ) 
f =i 

where Nj is the number of quark flavors and Sym is the Euclidean Yang-Mills 
action. The Euclidean four-volume is V4 = V3/?. The are non-Abelian gauge 
fields, which can be represented as 



T 

A - A a — 

I 1 \i 2 ' 



(6) 
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where the T a are the generators of the gauge group SU(iV c ). The number of colors 
is denoted by N c . The QCD Dirac operator is given by 

V = lfl (d tl + iA fl ). (7) 

This operator is anti-Hermitian, Tft = —T>. The 7^ are Euclidean gamma ma- 
trices with {7^)7v} = 23[j,v We use the chiral representation in which 75 = 
71727374 = diag (1, 1, -1, -1). 



2.2 Global Symmetries 

The structure of the QCD Lagrangian is to a large extent determined by sym- 
metries and renormalizability. As noted above, it is important to analyze these 
symmetries to construct the correct effective low-energy theory and the correct 
random matrix model. We now discuss three important global symmetries of the 
partition function and the Dirac operator. 



2.2.1 CHIRAL SYMMETRY AND TOPOLOGY 
The Dirac operator satisfies 

{75,^} = 0. (8) 

This relation is a compact expression of chiral symmetry, i.e. of the fact that 
right-handed and left-handed quarks can be rotated independently. One can 
write down an eigenvalue equation for T>, 

Vip n = i\ n ip n , (9) 

where the eigenvalues and eigenfunctions depend on the gauge field in Eq. (0). 
Using Eq. (||) one can show that the nonzero eigenvalues of T> occur in pairs ±iA n 
with eigenfunctions tp n and 75 ^n- There can also be eigenvalues equal to zero, 
A n = 0. The corresponding eigenfunctions can be arranged to be simultaneous 
eigenfunctions of 75 with eigenvalue ±1, i.e. these states have definite chirality. 
Denoting the number of zero eigenvalues with positive and negative chirality 
by N+ and iV_, respectively, the Atiyah-Singer index theorem states that v 
Y . Y is a topological invariant that does not change under continuous changes 
of the gauge field. However, the individual numbers N + and iV_ are not protected 
by topology, i.e. very small deformations of the gauge field will lift accidental zero 
modes. Thus, unless we impose very special constraints on the gauge fields, we 
generically have either N+ = or AL. = 0. 

In a chiral basis with r ysip R / L = ±ip R / L , one can use Eq. (Mj to show that 
( , ip R l \'D\ipn) = = ('i/'^J^IV'n) f° r an m an d n, where ^ = V^7o- From this 
property and the fact that T> is anti-Hermitian, it follows that the Dirac operator 
has the matrix structure 

(«•) 

This off-diagonal block structure is characteristic for systems with chiral symme- 
try. If there are n + v right-handed and n left-handed modes, the matrix W has 
dimension (n + u) x n, and the matrix T> in Eq. ( |i~0|) has \v\ eigenvalues equal to 
zero. 
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The QCD partition function can be decomposed into sectors of definite topo- 
logical charge u, 

oo 

Z QCB {0)= g e^QCD. 



v=— oo 



Here, the 0-angle is the coefficient of the topological FF term (which violates P 
and CP conservation) in the Lagrangian. We now show that the ^-dependence 
of the partition function can be absorbed into the phase of the quark masses. To 
this end, we introduce complex masses my for the right-handed quarks and the 
complex conjugate masses m*j- for the left-handed quarks. The partition function 
then reads 

^ QCD = (n^/' E(^ + m f m})) u , (12) 

/ A n >0 

where rhf = mj (m* f ) for v > (y < 0). In Eq. @), the average is only over the 
gauge field configurations with topological charge v, weighted by the Yang-Mills 
action. Since 



>*n™/'H r^T'Vi- 7 ( 13 ) 



f 



(e ie T\ f mA " for v > 



(e ie X\fm^j for v < , 



it follows that the ^-dependence of Z^ CD is entirely determined by the combina- 
tion e lS ]1/ m f ■ 

2.2.2 FLAVOR SYMMETRIES 

The second global symmetry is flavor symmetry, whose spontaneous breaking has 
profound implications for the hadron spectrum. To make this symmetry explicit, 
we rewrite the fermion determinant in Z^ CD as a Grassmann integral, 



N f r N f 

Y[ det(V + m f ) = dipdrp exp d' l xJ2 7 Pf(' D + m f) 7 P.f > ( 14 ) 

/=i L i=i 



where ip = V^7o- Again going to a chiral basis with r y^ R ' L = ±t/j R ' L , the 
exponent on the right-hand side of Eq. (FL4|) can be rewritten as 



r f - 

/ d 4 x J2 (ipfiW^f + tfiWtf + ipfmfipf + ipfmfipf) . (15) 
/ i 

In the chiral limit where all mf = 0, the fermion determinant is invariant under 
the transformations 

The only condition on U and V is that their inverses must exist. If the number 
of right-handed states Nr is equal to the number of left-handed states Nl, the 
symmetry is thus[| U(iV/) x \J(Nf). However, if Nr / Nl, i.e. if v ^ 0, the 

1 Actually, the symmetry group is Gl(iV/) x G1(7V/). The restriction to XJ(Nf) is a consequence 
of the Riemannian nature of the integration manifold (see Sec. ft.Sj). 
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axial symmetry group Ua(1) is broken explicitly by instantons or the anomaly. 
A second U(l) group, Uy(l), corresponds to the conservation of baryon number. 
The full flavor symmetry group in the chiral limit is thus given by G = SU(iVj) x 
SU(iV». 

If the quark masses are nonzero, the axial SU(iVy) subgroup with U = V^ 1 is 
broken explicitly by the mass term. The SU (Nf) vector symmetry (with U = V) 
is good for degenerate quark masses (mj = m for all /) but is broken explicitly 
for different quark masses {rut ^ my for all / 7^ /'). 

What is much more important than the explicit breaking, however, is the spon- 
taneous breaking of the axial-flavor symmetry. For an axial-flavor-symmetric 
ground state, the vacuum expectation value ('tpip) = (V^V^) + (V' i V ,R ) would be 
zero. However, phenomenological arguments and lattice QCD simulations indi- 
cate that ("tpi/j) ~ — (240 MeV) 3 . The spontaneous breaking of the axial symmetry 
also follows from the absence of parity doublets and the presence of Goldstone 
bosons, the pions. The quantity {ipip) is only invariant if U = V, i.e. the vac- 
uum state is symmetric under the flavor group H = SVy(Nf). Thus, the vector 
symmetries are unbroken whereas the axial symmetries are maximally broken. 
The first of these statements is in agreement with the Vafa-Witten theorem j2l| , 
which states that in vector-like theories such as QCD, vector symmetries cannot 
be spontaneously broken. The reasons why the axial symmetries are maximally 
broken [22] are less well understood. The Goldstone manifold is G/H = SU(iVy), 
and so there are Nj — 1 Goldstone bosons. 



2.2.3 ANTI-UNITARY SYMMETRIES 

Third, we consider the anti-unitary symmetries of the Dirac operator. Accord- 
ing to a fundamental theorem by Wigner, a symmetry in quantum mechanics is 
either unitary or anti-unitary. An anti-unitary symmetry operator A can always 
be written as A = UK, where U is a unitary operator and K is the complex 
conjugation operator. Below we always consider spectra of an irreducible sub- 
space of the unitary symmetries. If A = UK is an anti- unitary symmetry of the 
Dirac operator, then the symmetry operator A 2 = {UK) 2 = UU* is unitary, and 
in an irreducible subspace it is necessarily a multiple of the identity, UU* = Al. 
Because of this relation, U and U* commute so that A is real. By unitarity we 
have |A| = 1, which yields A = ±1. Therefore, the anti-unitary symmetries can be 
classified according to the sign of A 2 . There are three possibilities for the classifi- 
cation of Hermitian (or anti-Hermitian) operators: (a) there are no anti-unitary 
symmetries (denoted by the Dyson index (3 = 2); (b) if A 2 = 1, it is possible 
to construct a basis in which the operator is real (denoted by the Dyson index 
(3 = 1); (c) if A 2 = —1, it is possible to construct a basis in which the matrix 
elements of the operator can be organized into real (or self-dual) quaternions 
(denoted by the Dyson index /3 = 4). 

For QCD with iV c > 3 and fermions in the fundamental representation of the 
gauge group, there are no anti- unitary operators that commute with the Dirac 
operator. This means that the matrix W in Eq. ( |i~0|) is a general complex matrix 
with no further symmetries. There are two cases with nontrivial anti-unitary 
symmetries: QCD with two colors and fermions in the fundamental representa- 
tion, and QCD with any number of colors and adjoint fermions. We will discuss 
them next. 

For QCD with two colors and fermions in the fundamental representation, the 
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Dirac operator is given by Eq. (^) with = A a ^r a /2, where the r a are the SU(2) 
Pauli matrices in color space. The anti-unitary symmetry in this case is |23], ^4|] 

[Ct 2 K,W]=0, (17) 

where C = 7274 is the charge conjugation matrix. The square of the anti-unitary 
operator is (CT2K) 2 = 1. In this case, it is possible to find a basis in which the 
matrix W in Eq. ( |i~0| ) is real for all gauge field configurations [{L6||. 

As a consequence of the pseudoreality of SU(2), the symmetry of the QCD 



partition function in the chiral limit is enlarged to SU(2JV/) |22| , 25, ^6[. An 
axial U(l) is broken explicitly by instantons or the axial anomaly, as for N c > 3. 
The chiral condensate is only invariant under an Sp(2Nf) subgroup of SU(2iVy). 
The Vafa-Witten theorem, prohibiting the spontaneous breaking of global vector 
symmetries and assuming maximum breaking of the axial symmetries, thus pre- 
dicts a pattern of spontaneous chiral symmetry breaking given by SU(2iVj) — ► 
Sp(2iV / ). The Goldstone manifold is the coset SU(2iV / )/Sp(2A r / ), which is the 
set of antisymmetric unitary matrices. We thus have 2N 2 — Nf — 1 Goldstone 



bosons M, 26]. 



For fermions in the adjoint representation of the gauge group and any number 
of colors, the Dirac operator is = ^^(8^5^ + fabc-AV), where a and b are 
color indices and the f a b c are the structure constants of SU(iV c ). In this case, the 
anti-unitary symmetry is [^, |24| 

[CK,W]=0. (18) 

Because of (CK) 2 = —1, it follows that the eigenvalues of the Dirac operator 
are twofold degenerate with linearly independent eigenfunctions. In this case, 
it is possible to choose a basis in which the matrix elements of W are real (or 
self-dual) quaternions for all gauge fields. The eigenvalues of such a matrix are 
unit quaternions and therefore doubly degenerate |27]] in the representation of 
the Dirac matrix as complex numbers. 

Restricting ourselves to the case of even Nf, we can show that for Nf Majorana 
fermions, the flavor symmetry group is now SU(iVy). In this case, the chiral 
condensate is only invariant under an O(Nj) subgroup of SU(Afy). Applying 
again the Vafa-Witten theorem with maximum breaking of axial symmetry, we 
expect the pattern of spontaneous chiral symmetry breaking according to SV(Nf) 
— > 0(Nf), with the Goldstone manifold given by the coset S\J(Nf)/0(Nf). This 
is the set of symmetric unitary matrices. We thus have (Nf + 2)(Nf — l)/2 



Goldstone bosons [25, 26 1 



On the lattice, the symmetries of the Dirac operator may be different from the 



continuum symmetries [28]. We return to this point in Sec. 5.2. 



2.3 Dirac Spectrum 

Based on the eigenvalue equation (§) , we define the spectral density of the Dirac 
operator by 

P(A) = (£<H A - A ")}> (19) 

n 

where the average is over gauge fields weighted by the full QCD action. The 
spectral density is important because of its relation with the order parameter for 
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spontaneous chiral symmetry breaking, the chiral condensate {t^ip). It was shown 
by Banks and Casher [29] that 



E = |W>| = ^). (20) 

To be precise, we should have written S = lim e ^o hm m ^o limy^oo "xpi^/V, 
where it is important that the limits are taken in the order indicated. (In the 
normalization of Eq. (|l~9| ) the spectral density is proportional to the volume, so 
the explicit factor of 1/V in (f20|) is canceled to yield a finite result.) 

The relation (|20| ) can readily be derived. The chiral condensate is given by 

Id 

M) = - lim lim — — — log Z QCD (m) . (21) 
Nrr/ m^o v->oo VN f dm y ; v ' 

From Eq. (|5|), this yields 

(^M = - lim lim i\- r Y — — ) . (22) 

Since the nonzero eigenvalues occur in pairs ±iA n , their contribution to the sum 
can be written as 2m/(A 2 + m 2 ) (with \ n > 0). For gauge fields with topological 
charge u, the zero modes contribute a term |i/|/(mV). Assuming (u 2 ) oc V, 
we can drop these contributions in the limit V — > oo. In the same limit, the 
sum in Eq. ([22]) can be converted to an integral. In the limit m — > 0, we have 



2m/(A 2 + m 2 ) — > vr<5(A), which yields Eq. (|20|). As discussed above in Sec. LI 
spontaneous chiral symmetry breaking is encoded in an accumulation of the small 
Dirac eigenvalues; for the order parameter to be nonzero, we need p(0)/V > 0. 

An immediate consequence of the Banks-Casher relation is that the small eigen- 
values are spaced as 

AA = ^ = T^ < 23 > 

provided that p(0)/V > 0. This naturally defines a scale 

z = AT/S (24) 

for the study of the distribution of individual eigenvalues. For this purpose, it is 
convenient to define the so-called microscopic spectral density [20] 



*W = v5S.v^'(i^)- (25) 

This function describes the extreme infrared properties of the Dirac spectrum. 
Based on the arguments in the introduction, we expect it to be completely deter- 
mined by the global symmetries of the Dirac operator. As we demonstrate below, 
p s {z) can be computed both from the low-energy effective theory and from chiral 
RMT; the results coincide. Further confirmation of these ideas will come from 
results of lattice QCD simulations. 

At this point we would like to add some remarks on possible ultraviolet diver- 
gences. As an example, we consider the chiral condensate in the V — > oo limit, 
but before the limit m — > is taken. Converting the sum in Eq. ( p2|) to an 
integral and dropping the contribution from the zero eigenvalues, we obtain 

(#) = - lim lim [ A d\ 2 ™ P{X), V , (26) 
V Yl A^oom-^OJo A 2 +m 2 
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where we have introduced an ultraviolet cutoff A that must be removed at the end 
of the calculation. Asymptotically, the spectral density behaves as p(X) ~ VX 3 , 
which means that the integral is ultraviolet-divergent. However, this divergence 
does not contribute to {ipip} if the m — > limit is taken first to yield the usual 
Banks-Casher relation. Alternatively, we can first subtract the divergent contri- 



butions to the integral (26) and remove the cutoff before taking the chiral limit. 

A slightly different situation arises if we consider the dependence of the chiral 
condensate on a valence (or spectral) mass m v that does not appear in the average, 
so that p{X) is independent of m v . The quantity 

ES)= lim / Vy, (A)/ 7 (27) 



A^oo Jo X 2 + m 2 

is subject to ultraviolet divergences. Using again the fact that p(X) ~ VX 3 for 
large A, we find that the leading divergence is ~ m v A 2 . If we consider valence 
masses on the microscopic scale ( |24| ) we have m v ~ 1/(VS), and if the V — ► oo 
limit is taken before the A — > oo limit, the ultraviolet divergences are removed. 

In lattice QCD simulations with lattice spacing a, the cutoff is A = 1 /a. In the 
continuum limit, a — > 0, both the coupling constant and the ipvp-oper&tor have to 
be renormalized. For a discussion of the ultraviolet divergences in a recent lattice 
study of E(m v ) at finite volume, finite quark mass, and finite lattice spacing, see 
Ref. |0|. 



3 CHIRAL RANDOM MATRIX THEORY 

3. 1 Introduction of the Model 

In this section we introduce a chiral random matrix theory (chRMT) with the 
global symmetries of the QCD Dirac operator. In the spirit of the invariant 
random matrix ensembles, we construct a model with eigenfunctions distributed 
uniformly over the unitary unit sphere. This is achieved by choosing Gaussian- 
distributed random matrix elements. We thus arrive at the following chRMT 

Nf 

Z* (mi,...,m Nf )= / wndet(P + m/ )e-^ tff >, (28) 
where (5 is the Dyson index, 



, iWt o ) ' 

and W is an n x m matrix with v = m — n and N = n + m. The interpretation 
of this model is that N low-lying modes interact via a random interaction. A 
natural representation of this model is in the form of gauge field configurations 
given by a liquid of instantons. Then the low- lying modes are the zero modes 
of each instanton. We assume that v does not exceed \f~N so that, to a good 
approximation, n = N/2 for large N. The parameter N is identified as the di- 
mensionless volume of spacetime. For the formulation of this model with explicit 
factors of N/V included, see Ref. ||31[J . The potential v is defined by 



t#) = a ^ k ■ ( 3 °) 



k>l 
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The simplest case is the Gaussian case, where v{(f) = T, 2 (f). It can be shown (see 
Sec. |3.5| ) that the microscopic spectral density does not depend on the higher- 
order terms in this potential provided that the average spectral density near zero 
remains nonzero. The matrix elements of W are either real \fl = 1, chiral Gaus- 
sian Orthogonal Ensemble (chGOE)], complex [(3 = 2, chiral Gaussian Unitary 
Ensemble (chGUE)], or quaternion real \fl = 4, chiral Gaussian Symplectic En- 
semble (chGSE)]. In the latter case, the eigenvalues of V are doubly degenerate, 
and the use of Majorana fermions is implemented by replacing the determinant 
by its square root. For a non-Gaussian potential v(4>), we will omit the G in the 
abbreviations and use chOE, chUE, and chSE, respectively. Two earlier attempts 
to describe QCD Dirac eigenvalues used the Wigner-Dyson ensembles instead of 



the above chiral ensembles [32] 



This model reproduces the following symmetries of the QCD partition function: 

• The Ua(1) symmetry. All eigenvalues of the random matrix Dirac operator 
occur in pairs ±i\ n or are zero. 

• The topological structure of the QCD partition function. The Dirac matrix 
has exactly \v\ = \n — m\ zero eigenvalues. This identifies v as the topological 
sector of the model. 

• The flavor symmetry, which is the same as in QCD. For = 2 it is SU(iVj) x 
SU(JV», for (3 = 1 it is SU(2iV/), and for (3 = 4 it is SU(iV» (each Majorana 
flavor counts as 1/2 Dirac flavor). 

• The chiral symmetry, which is broken spontaneously with a chiral condensate 
given by 

S = lim tt P (0)/N . (31) 

N^oo 

(N is interpreted as the dimensionless volume of spacetime.) The symmetry- 
breaking pattern is g| SU(JV» x SU(JV» -» SU(JV», SU(2JV» -» Sp(2iV», 
and SU(A r j) — * O(Nf) for = 2, 1, and 4, respectively — the same as in QCD 



• The anti- unitary symmetries. These are implemented by choosing the matrix 
elements of W to be real, complex, or quaternion real for (3 = 1, (3 = 2, and 
f3 = 4, respectively. 

Along with the invariant random matrix ensembles, the chiral ensembles are part 
of a larger classification scheme that also includes ensembles for the description of 
disordered superconductors |3j|. In total, 10 different families of random matrix 
ensembles have been identified. They correspond one-to-one with the Cartan 
classification of symmetric spaces |3jJ]. 

The uniform distribution of the eigenfunctions over the unitary unit sphere is 
expressed as the invariance 

W -► rfWV , (32) 

where the n x n matrix U and the m x m matrix V are orthogonal matrices 
for (3 = 1, unitary matrices for (3 = 2, and symplectic matrices for (3 = 4. 
This invariance makes it possible to express the partition function in terms of 
eigenvalues of W defined by 

W = U ] KV . (33) 
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Here, A is a diagonal matrix with real diagonal matrix elements > 0. In terms 



of the eigenvalues, the partition function (28) is given by 



Z^( mi , . . . , mNf ) = [ d\\A(X*)f J] 4 V^2> nJfiXl + rn}), 

f 

(34) 



k 

where the Vandermonde determinant is defined by 



A(A 2 ) = n(Al"Af). (35) 
k<i 

In Eq. ( p4[) and elsewhere in this review, we have omitted the normalization 
constant of the partition function. From the joint eigenvalue distribution (which 
is the integrand of Eq. (|34l) ), we see that a nonzero topological charge v can be 
introduced by adding j3\v\/2 massless flavors to the theory with v = 0. Therefore, 
this duality between flavor and topology is a general feature of all correlation 
functions. For a discussion of this duality in terms of finite volume partition 



functions, see Refs. |35, 3q| . 



3.2 Sigma Model Representation of Chiral Random Matrix Models 

The chiral random matrix partition function can be evaluated using standard 
random matrix methods. For simplicity, let us consider the case (3 = 2. The 
fermion determinant can be written as a Grassmann integral, and averaging over 
the Gaussian distribution function results in a four- fermion interaction. Because 



of the underlying unitary invariance (32) of the chRMT, the fermionic variables 
only appear in invariant combinations of the form o^ g ~ Wipf, and the partition 



function can be rewritten identically in terms of these variables as [p0| , |3?j| 

Z u {m)= J Dae~ nS2TTaat det u (a + ml)det n [(a + ml)(a^ + ml)\ , (36) 

where the quark masses have been taken real and degenerate for simplicity. The 
integration is over the real and imaginary parts of the Nf x Nf arbitrary complex 
matrix a. For m = this integral is invariant under a — ► UaV^ 1 with U and 
V £ SU(A^j). Let us calculate the integrals by a saddle-point approximation. 
The saddle-point equation reads (notice that v -C n) 

S 2 (a f + m)<r = 1 . (37) 

The condensate at the saddle point is given by 

= 7^ W d m \ogZ= JL(Tr(<7 + at)) , (38) 



where the expectation value is with respect to the partition function (36). The 
solution of the saddle-point equation, a, is proportional to the identity matrix. 
In the chiral limit, we then find a = 1/S, and we can identify £ as the chiral 
condensate. 

The fluctuations about the saddle point can be separated into massive modes 
with a curvature of order n and Goldstone modes with a curvature of order ran, 
In the limit n — ► do, the partition function (|36|) can be simplified by keeping the 
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integrals over the Goldstone manifold and performing the remaining integrals by 
a saddle-point integration. This results in the partition function 

Z v (m) = [ DU det»U e nSTr i^+MW^) ; (3g) 

JU&J(Nf) 



which is the familiar finite volume partition function to be discussed in Sec. 4.1 



3.3 8 -Dependence of the Partition Function 

The ^-dependence of the partition function is obtained by summing over all topo- 
logical sectors of the partition function according to Eq. (|ll~|), which applies to 
full QCD where the partition function is obtained by integrating over the gauge 
fields. The terms Z v in this equation reflect the probability of encountering gauge 
fields with topological charge v. 

For the effective theory, the partition function in the topological sector v is 
given in Eq. (|3^). However, this ^-dependence is only due to the fermion deter- 
minant. In the sum over v, we should therefore add weight factors P{v) that take 
into account the distribution of topological charge of the quenched gauge fields 
(i.e. gauge fields generated without the fermion determinant). This yields 

Z eS (9, m) = J2 P{v)e iv9 Zf(m) . (40) 

V 

Below, we argue that the weight factors can be ignored for light quarks but that 
they are necessary when one considers the quenched theory (which corresponds 
to Nj = or, equivalently, to the limit of very heavy quarks). 

Our starting point is that the topological susceptibility of the quenched gauge 
fields is nonvanishing, i.e. 

(41) 

Invoking the central limit theorem, we assume that P{v) is Gaussian, 

1 



P(v) = . e 2 <^ , (42) 

which has been verified by quenched lattice QCD simulations. 

We analyze the partition function ( |40"|) for Z„(m) given by the finite- volume 
partition function in Eq. (|3^) . Replacing the sum over v by an integral, we obtain 
after integration 

reBfa m \ _ f nTT -■^ k (0-iTrlo g U) 2 +iEyTr(Mt/+Mta- 1 ) 



The exponent in the integrand corresponds to the effective static potential of QCD 
with a large number of colors pS[ ], and therefore the chRMT partition function 
reproduces all identities that have been derived from this effective potential ]20| , 
|39| ]. For example, the topological charge is screened by light quarks. For a 
careful analysis of the periodicity requirements of the ^-dependence of the effective 



partition function (43), see Ref. |4( 

Let us illustrate the screening of topological charge for Nf = 1 and real m, for 
which 



Z cS (6,m) = I #e-^( 9 -« 2+SVmco ^ . (44) 
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The integral over <j) can be performed by a saddle-point approximation resulting 

Z cS (d,m) = e mVJ]cose . (45) 



m 



The topological susceptibility is given by 

d 2 

(z, 2 ) = _ log Z = mVY, (for 6 = 0) (46) 

Oo 

so that the topological charge is completely screened in the chiral limit. This also 
shows that the results are insensitive to our choice of the distribution function 
P(v). For example, in an instanton liquid interpretation of the chiral random 
matrix model, the natural choice for P{y) is the binomial distribution B(N —v, v), 
resulting in the same effective potential p0| . 

In fact, it can be concluded from Eq. (|46|) that for light quarks the distribution 
P(v) can be ignored altogether. Let us again illustrate this for Nf = 1, for which 
the integral in Eq. ( |39| ) results in the finite volume partition function 

Zf(m) = I u (mVZ) , (47) 



where I u is a modified Bessel function. Summing over v according to Eq. (|T 
i.e. without the weight factors P{v), we obtain plf 

Z eS (9, m) = J2 &ivt } Iu{mVT) = e mVScos0 , (48) 

V 

where we have used a summation formula for modified Bessel functions. This is 
exactly the same result as in Eq. (|45|). For an extension of these results to more 
than one flavor, see Ref. [41]. 



3.4 Spectral Correlation Functions 

In this section, we discuss the spectral correlation functions and the microscopic 
spectral density corresponding to the chRMT partition function. The global 
spectral density is a semicircle. It is not universal and we do not expect to find 
it in realistic physical systems. The universal quantities are local or microscopic 
spectral correlation functions. In general, a fc-point spectral correlation function 



is defined by [27 



R k (X 1 ,...,X k ) = / Y, S ( X i- E h)--- S ( X k-E jk )\ , (49) 

jp ¥= jq 

where the Ej p are eigenvalues. The quantity R k represents the probability density 
to find k eigenvalues, regardless of labeling, at Ai, . . . , In particular, Ri(X) = 
p(A). By "local" we mean that the energy differences |A P — A 9 | are of the order 
of a few mean level spacings. 

Because of the Ua(1) symmetry of QCD, all nonzero eigenvalues come in pairs 
±iX n (see Sec. |2.2.1| ). This implies that the origin, A = 0, is a special point 
of the spectrum and that we have to consider the R k separately in the bulk of 
the spectrum and near A = 0. The latter is called the microscopic region or the 
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"hard edge" of the spectrum.^] In addition, there is the tail or "soft edge" of 
the spectrum, which we do not discuss. The bulk of the spectrum is the middle 
region, far from either edge. 

Because the random matrix partition function is known in terms of an integral 
over the eigenvalues of the Dirac operator, see Eq. it is possible to obtain 

the spectral density and all spectral fc-point correlation functions by integration 
over N — k eigenvalues. In the present case, the integrals can be performed most 
conveniently by the orthogonal polynomial method. We mention only the most 
important results in the sections below. 



3.4.1 BULK CORRELATIONS 

Our conventions are such that if we have N eigenvalues, the global spectral 
density is normalized to N so that the mean level spacing is of order 1/N. Local 
spectral correlation functions R k (Xi, • • • , A/%) in the bulk of the spectrum are thus 
characterized by A p ~ 0(1) and |A P — A J ~ 0(1/N). Universal correlation 
functions are obtained by rescaling the eigenvalues according to the average local 
level spacing. This procedure, known as unfolding, is discussed in more detail in 



Sec. |5Tj. The rescaled eigenvalues with average level spacing equal to unity are 
denoted by x k - 

In the bulk of the spectrum, the RMT results for the R k are identical for 
the chiral ensembles and the corresponding nonchiral ensembles For 
simplicity, let us consider the simplest ensemble, the chUE. On the unfolded 
scale, we obtain in the limit N — > oo 

R k (xi,. . . ,x k ) = det[K(xp,x q )] Ptq=1) ... }k (50) 

with the sine kernel 

smir(x-y) 

K{x, y) = - f r- • (51 

vr(x - y) 

The fact that K(x,y) depends only on \x — y\ reflects the translational invariance 
of the spectral properties after unfolding. The results for the chOE and the chSE 



are somewhat more complicated than those for the chUE [17]. 

Other quantities can be derived from the R k . As a measure of short-range 
correlations between the eigenvalues, one considers the nearest-neighbor spacing 
distribution P(s), which is the probability density to have a spacing of s be- 
tween adjacent levels. In random matrix theories, one finds that P(s) ~ for 
small values of s and has a Gaussian tail for large spacings. Long-range spectral 
correlations are characterized by the level- number variance S 2 (L) given by 

£ 2 (L) = ((n(L)-L) 2 ) , (52) 

where n{L) is the number of levels in an interval of length L, or by the spectral 
rigidity A^L) f44)] , which is defined as an integral transform of S 2 (L), 



2 



L 



A 3 (L) = -j du(L-u)(L 2 - Lu-u 2 )Z 2 (u) . (53) 
Lr Jo 



2 The reason for this nomenclature is that in a Fermi-gas formulation of RMT, where the 
eigenvalues are interpreted as the positions of particles, this symmetry corresponds to choosing 
a potential that is infinite for negative values of the eigenvalues. In this picture, the soft edge 
corresponds to a potential with a finite slope. 
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The advantage of using the A3-statistic is that its statistical fluctuations are 
much smaller than those of the £ 2 -statistic. The quantities X 2 and A3 can be 
derived from the two-point function R2. To compute P(s), one needs all fc-point 
functions. These quantities are very different for systems with correlated eigen- 
values (described by RMT) and systems with uncorrelated eigenvalues (described 
by the Poisson ensemble). For example, the large- L behavior of the chUE results 
is given by S 2 (L) ~ (lnL)/-7r 2 and As(L) ~ (lnL)/27r 2 , whereas for uncorrelated 
eigenvalues one finds X 2 (L) = L and A%(L) = L/15. 



3.4.2 MICROSCOPIC CORRELATIONS 

In the context of QCD, the small eigenvalues are more interesting than those in 
the bulk because of their relation to spontaneous chiral symmetry breaking via 
the Banks-Casher relation. The local spectral correlation functions at the hard 
edge of the spectrum are characterized by X p ~ 0(1/N) and \X P — X q \ ~ 0{1/N). 
Unfolding in the microscopic region is easy; the eigenvalues are simply rescaled 
by the mean level spacing at A = 0, which is given by Eq. (p3|). By convention, 
we define the unfolded variables by x p = X p VT,, i.e. the factor of n is omitted. 
Note that V oc N. Thus, in the unfolded variables we have x p ~ 0(1) and 
\xp — x q \ ~ 0(1). As advertised above, the functional form of the Rk is different 
in the microscopic region. The Rk are still given by the determinant of a kernel 
according to Eq. (|50|), but the sine kernel is now replaced by the Bessel kernel 



[45] 



x z — y z 

where J a denotes Bessel functions and a = Nf + \ v\, with Nf the number of 
massless flavors and v the topological charge. A mathematical discussion of the 



Bessel kernel can be found in Ref. [46] 



The RMT results for the chOE and the chSE are more complicated 



49, 33]. Note that, in contrast to the situation in the bulk, the microscopic 
correlations depend on Nf and v. 

The one-point function is universal in the microscopic region. For the chUE, 
the microscopic spectral density, defined above in Eq. (p5[), is related to the kernel 



by p s (z) = K(z, z) and can be expressed as |45, p0| 



Ps{z) = ^[Jn }Mv \{z) - J Nf+ \ v \ +1 {z)J Nf+ \ v \_ 1 {z)\ [ . (55) 

It is also the generating function for the Leutwyler-Smilga sum rules, which are 
essentially inverse moments of p s (or of higher-order spectral correlation func- 
tions). They are discussed below. The results (^4j) and ( pq ) were first derived 
using chiral RMT. In the meantime, Eq. ( |55| ) has been confirmed by an explicit 
calculation starting from a partially quenched chiral Lagrangian |3l], p2; ], which is 
discussed in the next section. The microscopic spectral density and the two-point 
correlation function for the case Nf + v\ = have also been derived by means of 
the supersymmetric method of RMT 53]. 

There is a very interesting connection between the Bessel kernel and the ef- 
fective partition function. In addition to the partition function Z^ N f^ for Nf 
massless sea quarks, one can also compute the partition function Z^ N f +2 ^ for Nf 
massless quarks and two additional flavors with imaginary masses. One then 
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finds the relation 54] 

K(x,y) = -(xy) i + 1 , (56) 

which can also be generalized to massive sea quarks. Thus, to obtain the Rk 
for Nf flavors, one needs to know the partition function for two additional fla- 
vors with imaginary masses. This result has a natural interpretation in terms 
of the effective spectral partition function (discussed below in Sec. |). The two 
additional flavors are the spectral quark and its bosonic superpartner The 
relation ( p6| ) once again demonstrates the equivalence of the random matrix for- 
mulation and the effective field theory in the zero mode domain. Additional 
consistency conditions for the finite volume partition function are discussed in 
Ref. 



3.5 Universality Proofs 



In most random matrix calculations, the probability distribution of the random 
matrix elements is assumed to be Gaussian to simplify the calculation. Nonuni- 
versal features such as the global spectral density depend on the choice of the 
probability distribution. On the other hand, universal results of RMT do not 
depend on this choice, nor on other deformations of the random matrix model. 
In addition to numerous empirical verifications of universal behavior in RMT, 
analytical proofs of universality have recently been constructed for two different 
types of deformations. 

First, consider deformations that do not break the invariance of the theory 
under unitary transformations of the random matrix. In this case, the Gaus- 
sian distribution P(W) oc exp(— Na\ Tr WW* 1 ) is replaced by the more general 
distribution P(W) oc exp(— iV J2T=i a k Tr (W W>) ), i.e. the quadratic term in 
the exponent is replaced by an arbitrary polynomial. The advantage of the uni- 
tary invariance is that the probability distribution depends only on the eigen- 
values of WW^ . For the chiral ensembles, it has been shown that both the 
bulk and the microscopic spectral correlations remain unchanged on the un- 
folded scale @, H, H H, HE H- The weight function of QCD also con- 
tains the fermion determinant, and it has been shown that the universal results 
are unchanged if the random Dirac matrix in the determinant is replaced by 
a polynomial of that matrix [61]. Universal microscopic correlations have also 
been found for the so-called Ginsparg- Wilson Dirac operator 62]. Along with 
work on universality for the much more widely studied Wigner-Dyson ensem- 

these findings make it clear 



bles Wc 



65, 66, 



6S, 57, WO, 71 



m 



that spectral correlations on the scale of the average level spacing are strongly 
universal, i.e. they do not change despite substantial variations of the average 
spectral density. There have also been several interesting results on the validity 
of universality for wide correlators, i.e. on macroscopic scales | 



Second, consider deformations that do break the unitary invariance by adding 
an arbitrary deterministic matrix Y to the random matrix W in the Dirac oper- 
ator. As we discuss in Sec. |(| such a model can provide a schematic description 
of QCD at nonzero temperature. In particular, the matrix Y can be chosen in 
such a way that p(0)/V vanishes so that chiral symmetry is restored. It has been 
proven that the bulk correlations on the unfolded scale are unaffected by the 
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matrix Y, and that the same is true for the microscopic spectral correlations as 
long as p(0)/V remains nonzero pi WK W% [7|, |7|, 



Most universality proofs have been performed for random matrix ensembles 
with (3 = 2. However, it is possible to establish relations between the (3 = 
2 kernel and the kernels of the other two ensembles with (3 = 1 and (3 = 4 
[47, 81, 82, 83, 84|. Prom these relations and the universality of the (3 = 2 results, 
one can then infer the universality of the (3 = 1 and (3 = 4 results. Recently, 
interesting connections between massive correlators and the massless correlators 
for (3 = 1 and (3 = 4 have been made [85], and the transition between these 
ensembles and the (3 = 2 ensemble has been studied [86]. 

The microscopic spectral correlations are universally given by RMT only if 
p(0)/V > 0. An interesting question [^7]] is what happens to the Rk if there is 
some phase transition so that p(0)/V vanishes. This transition could be triggered 
by, for example, nonzero temperature or a large number of flavors. Beyond the 
transition point, there will be a gap in the spectrum, and the smallest eigenvalues 
beyond this gap are presumably described by the soft-edge results of RMT. What 
is more interesting, however, is the transition point at which p(0)/V becomes 
zero so that there is no gap. This situation has been investigated in two different 
ways. First, keeping the unitary invariance, one can fine-tune the polynomial in 
the exponent of P(W) to make p(0)/V just vanish [87]. Second, breaking the 
unitary invariance, one can choose a "critical" deterministic matrix Y so that 
p(0)/V 



without a gap 



In both cases, one obtains new functional 
forms for the microscopic spectral correlations and a different scaling with the 
volume (or matrix dimension). However, the results are not unique; they depend 
on the way one chooses to approach the p(0)/V = limit. It is, therefore, an open 
question as to whether the microscopic spectral correlations of the QCD Dirac 
operator at the chiral phase transition agree with one of the results obtained in 
a random matrix model. 



4 EFFECTIVE THEORIES AT LOW ENERGIES 

4-1 Finite- Volume Partition Function 

As we have seen before, chiral symmetry is spontaneously broken by the QCD 
vacuum. According to Goldstone's theorem, this leads to the appearance of 
massless modes. The Lagrangian describing the dynamics of these modes can be 
obtained solely from the symmetries of the QCD partition function. Goldstone 
modes corresponding to the symmetry-breaking pattern SUn(Nf) x SU^Afy) — > 
SUy(iVy) can be parameterized as 

U = UrUl 1 (57) 

with Ur € SUij(iVj) and Ul £ SUl(A^). The chiral Lagrangian is constructed 
by the requirement that it should have the same invariance properties as the 
QCD Lagrangian. In particular, for m = the Lagrangian should be both 
Lorentz invariant and invariant under SU r(Nj) x SU i{Nf). To lowest order 
in the momenta (or derivatives), the kinetic term in the effective Lagrangian is 
uniquely given by 0] 

Adn = ^TV d^Ud^U- 1 . (58) 
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The parameter F is the pion decay constant. The mass term in the QCD La- 
grangian breaks the full flavor symmetry. However, the full symmetry can be 
restored if the mass matrix is transformed as well, 

M UlMUk 1 . (59) 

We require that the mass term in the effective Lagrangian satisfies this extended 
symmetry. To lowest order in M, it is therefore uniquely given by 

C m = --£ Tr (MU + M^U~ X ) . (60) 

For a diagonal mass matrix, the action is therefore minimized by U = 1. The 
normalization of the mass term is such that the mass derivative of the partition 
function according to Eq. (21) is equal to the chiral condensate S. The chiral 
Lagrangian is valid in the domain 

m < A and p«A, (61) 

where p is the momentum and A is a typical hadronic mass scale. 

We will study the effective Lagrangian at finite volume in a box of volume L 4 . 
Then the smallest nonzero-momentum modes are of the order 

p~l/L. (62) 

The fluctuations of the zero-momentum modes, the constant fields, are not af- 
fected by the kinetic term and are limited only by the mass term. Comparing 
Eqs. fl58D and (p0|), we see that in a domain where 

mT, 1 

-pT « Jl > (63) 

the fluctuations of the zero-momentum modes completely dominate the fluctu- 
ations of the nonzero-momentum modes. The mass-dependence of the effective 
partition function is then given by [19, 23| 

Z(m,0)= f DU e VXReTrUme W/N f ^ ^ 

Juesv(Nf) 

where we have set M = ml and have introduced the ^-dependence by the sub- 
stitution m — > mexp(i8/Nf) just as for the QCD partition function. Combining 
the conditions ( |61~D and (|63[ ), we find that the finite volume partition function is 
valid in the domain |l^, E3[ 

1 

m < — r < A . (65) 
AL 2 

The partition function for a given topological charge can be extracted from its 
^-dependence. Since 

Z(m,6) = Y j e iue Zu{m) , (66) 

V 

we obtain Z v by Fourier inversion. Thus, the finite volume partition function in 
the sector of topological charge v is 

Z„(m) = — / d6e- iv6 Z(m,9) = [ DU det"U e VSKeTrmU , (67) 

27T J JU£V{N f ) 

which coincides with Eq. (|39|). The partition function fl67D can be evaluated 
for arbitrary masses. The integrals were first calculated in the context of one- 
plaquette lattice QCD models [9C] but later rederived by means of Itzykson-Zuber 
integrals M, 62, M. 
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4-2 Leutwyler-Smilga Sum Rules 



For masses in the domain of Eq. (|6q), we can equate the mass dependence of the 
QCD partition function and the mass dependence of the effective partition func- 
tion. Equating the coefficients of the expansion in powers of the mass then gives 
us the so-called Leutwyler-Smilga sum rules [23|. The QCD partition function 
can be expanded as 

z «cD (m) = '"' IH ni-H+^)> = mM (1 + E + - 

where the prime indicates that the product is over nonzero eigenvalues only. The 
effective partition function can be expanded in a power series in m as well. From 
the U(l) part of the group integral, it is clear that Zf(m) oc m)^. Thus, 

Zf = m^(a + a 1 m 2 + ...), (69) 

where the coefficients en are obtained by calculating the group integrals. Let us 
consider two examples. The simplest example is Nf = 1, for which 

Zf = — [ dd e -iv6 e mVX COS 8 = Iv ( mVT ,) . (70) 



2vr 

The series expansion of the modified Bessel function is given by 

for v > 0. We also have I- V (x) = I u (x). Matching the normalizations of Eqs. ( |6"8| ) 
and ( |69| ) and equating the coefficients of the 0{m 2 ) terms, we find 

(72) 




4(H + 1) 



A second example is the case of v = and arbitrary Nf. In this case, we use the 
group integral 

J DUU^ul^^Su (73) 



to derive the result 

^\£^) = ^v (74) 



Since asymptotically for large A the spectral density is proportional to V\ 3 , the 
sum over the eigenvalues has to be regularized. As discussed above in Sec. a 
finite result is obtained by taking the thermodynamic limit before removing the 
cutoff. 

What do we learn from these sum rules? Since the total number of eigenvalues 
is of order V, the obvious interpretation is that the smallest eigenvalues are of 
order 1/V. Indeed, this is in agreement with the Banks-Casher formula. What 
is more important, however, is that in the derivation of the sum rules we have 
relied only on the chiral symmetry of the partition function and its spontaneous 
breaking by the formation of a chiral condensate. Therefore any theory with 
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the same pattern of chiral symmetry breaking as QCD should obey the same 
spectral sum rules. In particular, the eigenvalues distributed according to the 
chRMT introduced in the previous section should obey the same sum rules. In 
fact, all Leutwyler-Smilga sum rules can be derived systematically from the joint 
eigenvalue chRMT probability distribution in Eq. ([34]) using the known Selberg 
integrals [94]. In addition to these sum rules, it is possible to derive massive sum 



rules |9f|. A relation between sum rules and partially quenched effective theories 
was considered in Ref. [P^]. Spectral sum rules for Ginsparg- Wilson fermions 



on the lattice were considered in Ref. [97]. Inverse moments of the small Dirac 



eigenvalues in generalized chiral perturbation theory were derived in Ref. 9Jj. 

We have seen that chRMT allows us to calculate the correlations of the smallest 
eigenvalues of the Dirac operator on the scale of the average level spacing. This 
raises the question of whether it is possible to derive the microscopic properties 
of the eigenvalues from the low-energy effective partition function. The answer 
is no. The spectral density cannot be derived from the mass dependence of the 
chiral condensate, 

Z(rn) = (±y:-^—-), (75) 



because the average contains the fermion determinant with the same mass. As we 
show in the next section, this problem can be circumvented by the introduction 
of a spectral mass. 

4-3 Spectral Mass and Partially Quenched Partition Function 
The Dirac spectrum can be obtained from the resolvent 

were the spectral mass z is an independent complex variable that does not occur 
in the average. With purely imaginary eigenvalues the spectral density is given 
by the discontinuity of the resolvent across the imaginary axis, 

^ = lini ii[S(iA + e) - E(*A - e)] . (77) 

This follows immediately upon writing the trace as a sum over the eigenvalues of 
T>. The generating function for the resolvent can be written down easily |99|| , 

_ / det w /(P + m)det(P + z) \ 

=\ ; ym - < 78 > 

where we have explicitly displayed the fermion determinant in the measure. The 
resolvent is then given by 



£(*)= ^d z Z^(z,z') 



(79) 

z 

In QCD, this method was first introduced in order to derive the quenched strong 



coupling expansion [10C[. In nuclear physics and condensed matter physics, this 



method is known as the supersymmetric method or Efetov method for quenched 



disorder [101, 102, 103]. 
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The determinants in the numerator can be written as fermionic integrals where- 
as the determinant in the denominator can be written as a bosonic integral. 
The partition function ( |7q ) is thus invariant under flavor symmetries that mix 
commuting and anticommuting degrees of freedom. More precisely, for m = z = 
z' = the partition function is invariant under the supergroup G\r{Nj + 1|1) x 
G\i{Nf + 1|1). In order to obtain a consistent effective theory, it is essential to 
extend the unitary symmetry to the general linear group Gl. We expect, as for 
the QCD vacuum, that chiral symmetry is spontaneously broken to the diagonal 
subgroup Q\y{Nf + 1|1). The mass term explicitly breaks the full symmetry to 
this subgroup as well. 

The effective partition function can be derived in the same way as the regular 
chiral Lagrangian with one complication. To obtain finite integrals, we must make 
sure that the integration manifold is Riemannian. This is why we have extended 
the flavor symmetry to the full general linear group. The integration manifold is 
then given by the maximum Riemannian submanifold of G1a(Aj + 1|1). In plain 
language, this means that we compensate the extra minus sign in the supertrace 
by complexifying the group parameters. The generating function is thus given 



by [51] 



Z pq = / DUe 



El- 

4 



Str d^Ud^U- 1 -^^ (-^-f) +f Str (MU+AtfiU- 1 ) 



U<EG\(N f +l\l) 

(80) 

where iV2&o/F = Strlogt/. This partition function is also known as the par- 
tially quenched effective partition function. It has been used to obtain a better 
understanding of quenched lattice QCD results [p. 04j, |105f| . The hat on Gl denotes 
the maximum Riemannian submanifold of Gl(Nf + 1|1), and Str stands for the 
supertrace. An example of an explicit parameterization of U will be given in the 
next subsection. The parameter mg is proportional to the topological suscepti- 
bility and results in a mass for the singlet channel that does not vanish in the 
chiral limit. 

The perturbative formulation of this partition function was first given by 
Bernard and Golterman 1 104 ] , and a perturbative calculation of the dependence 



of the condensate on the spectral mass was first performed in Ref. [ |106| ]. In that 
case, it is acceptable to ignore the convergence properties of the effective partition 
function and integrate over the unitary supergroup instead. 

An alternative generating function would be obtained by replacing the ratio of 
the two determinants by det n (P + z) and putting n — > after having calculated 



the resolvent. This procedure, known as the replica trick [107], successfully repro- 
duces the asymptotic expansions of the spectral correlation functions | |108] , |109| ]. 
There have been recent claims that it is possible to derive truly nonperturbative 



results by means of the replica trick 110, 111, 1 12 1 . We do not believe that 
this can be done in general glf, fiUf . As mentioned, because of the com- 



pact/noncompact structure of the final answer, the supersymmetric formulation 
is the natural approach to this problem. 

4-4 Domains of the Partially Quenched Effective Theory 
The Goldstone fields can be written as 

U = e l ^ u ' F . (81) 
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To second order in the fields, the effective Lagrangian in momentum space is 
given by 

£ = 77EE( fc2 + M a)n^), (82) 
a k 

where the sum is over the momenta in a box of length L (including the zero- 
momentum state) and the masses of the Goldstone bosons are denoted by M a . 
The sum over a is over the different Goldstone modes. In addition to the usual 
Goldstone modes with mass 

2 2mS 

M mm = p 2 ' \°^) 

there are both fermionic and bosonic Goldstone bosons with mass 

ML = (84) 
and fermionic and bosonic Goldstone bosons with mass 

Ml = 2 -§ ■ (85) 

In the Lagrangian (p2[), one can distinguish the zero- momentum modes from the 
nonzero-momentum modes. The magnitude of the fluctuations of the nonzero- 
momentum modes is of order 1 / (k 2 + M% ), whereas the magnitude of the fluctua- 
tions of the zero-momentum modes is of order 1 /M% . Since the smallest nonzero 
momenta are of order 1/L, the fluctuations of the zero- momentum modes are 
dominant if 

M 2 a « • (86) 

This means that the Compton wavelength of the "pion" is much larger than the 
size of the box. For nonzero quark masses of order 0(L°), the inequality (86) 



is never satisfied. However, the spectral mass z is a free parameter, and the 
inequality can be rewritten as [|116| , 117] 



z « E c = ^ . (87) 
We sometimes refer to the quantity E c as the Thouless energy for reasons that will 



become clear in Sec. 7.3. In the domain of Eq. (|87j), the dominant contributions to 
the resolvent are from the zero-momentum modes. Thus, the partially quenched 
effective partition function J8C| ) can be reduced to the partition function in the 
zero-momentum sector [51] ) |52| , 

Z= [ DUe hWStr (MU+MtU-i)_ (8g) 

JueGl(N f +l\l) 

The microscopic spectral density and the spectral correlations computed from 



Eq. (88) are identical to those obtained in chRMT |116| , 117 1 . Let us show this 
by an explicit calculation of the dependence of the chiral condensate on the 
spectral mass z. To cover the complete range (|87|), the integral over U has to 
be done nonperturbatively. This is a straightforward superintegral that can be 
evaluated using standard methods. The calculation of the integration measure 
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requires an explicit parameterization of the integration manifold. For example, 
in the quenched limit, Nf = 0, a possible choice is 



U 



e LV a 
e s 



with a and (3 Grassmann variables, <fi E [0, 2n], and s G (—00,00). A charac- 
teristic feature of the integration manifold is that it consists of a compact and 
a noncompact component. The valence-quark mass dependence of this partition 
function and its generalization to arbitrary topological charge coincide with the 
valence quark mass dependence of the chUE partition function. We merely quote 
the final result for the resolvent for Nf massless flavors in the sector of topological 
charge v 1116(1 , 

^ = u [I a (u)K a (u) + I a+l (u)K a ^{u)\ + M , (90) 

u 

where a = Nj + \v\ and u = zVH. The compact/noncompact symmetries are 
reflected in the appearance of the I a /K a -Bessel functions and are thus a natural 
ingredient of the underlying integration manifold. This result was first obtained 
from chRMT by integrating the microscopic spectral density in Eq. (^) as follows, 



E(w) _ f 00 2u 



<72-^2^(0 • (91) 



Alternatively, the microscopic spectral density can be obtained by taking the 
discontinuity of according to Eq. (fT7|). By integrating back, we find that 
the mass dependence of the zero- momentum partially quenched partition function 
coincides with that of the chRMT partition function for (3 = 2. 

Figure || presents a schematic picture of the different domains in the Dirac 
spectrum. There are other non-QCD theories that can be reduced to the same 
partially quenched partition function. Probably the best known example is the 
instanton liquid model of QCD p|. Another example, more closely related to 
disordered condensed matter systems, is the random-flux model |118|, which is 



in essence quenched lattice QCD with Kogut-Susskind fermions but without the 
phase factors due to the 7-matrices. Related examples are so-called two-sublattice 



models with disorder 119 ], and disordered lattice models with the chiral and 
flavor symmetries of QCD flj(], [121 . 
In the range 

E c < z <g A , (92) 

the effective action (p0|) is still valid but chiral random matrix theory no longer 
applies. In this range, the U fields can be expanded to second order in the pion 
fields, and the spectral mass dependence of the chiral condensate can be obtained 
to one loop order. The spectral density follows by taking the discontinuity of 
T,(z) across the imaginary axis. This calculation can be performed both for the 
partition function of Eq. (^(]) and the effective partition functions for (3=1 and 
(3 = 4. The result for the slope of the spectral density at zero, valid for Nf > 2, 



reads [|l_22 



P'(0) = {N f -2){N f + (3) S 
p(0) 16n(3N f F 4 ' K ' 

The result for (3 = 2 was first derived by Smilga and Stern |122j ] from the scalar 
susceptibility in standard chiral perturbation theory. We emphasize that this is 
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non-QCD 
chiral theories 



z < A 




v 



AL 2 



chRMT 



Figure 2: Diagram illustrating the different domains in the QCD Dirac spec- 
trum. Here, pqChPT denotes partially quenched chiral perturbation theory, and 
chRMT stands for chiral random matrix theory. 



an exact result for the QCD Dirac spectrum that is valid in the thermodynamic 
limit. The vanishing of the slope for Nt = 2 has been confirmed by instanton 



liquid simulations [123]. 

The two-point correlation function can also be obtained from a super symmetric 
generating function with, in this case, two different spectral masses and two 
bosonic superpartners [51|. Again the zero-momentum contribution coincides 
with the chUE result. Let us consider the perturbative expansion of the two- 
point correlation function. Generalizing Eq. (|77|), we obtain a relation between 



the two-point level correlation function and the pion susceptibility [51, 124], 



(p(A)p(A')) 



c 



4tt 2 F 4 



Disci 



z=i\,z'=i\ 



(q 2 + M 2 ) 2 ' 



(94) 



where the meson mass is given by M 2 = (v?+ \fzJ 2 )Y J jF 2 and the superscript 
C denotes the connected part of the two-point function. By taking the discon- 
tinuities of the zero-momentum contribution, we find the two-point correlation 



function 51] 



(A -A') 2 (A + A') 2 



(95) 
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which is the correct asymptotic result for the chUE. 
4-5 QCD in Three Euclidean Dimensions 

QCD in three Euclidean dimensions can be analyzed in much the same way as 
discussed above. The main difference is the absence of the Ua(1) symmetry and 
the absence of instantons. In terms of the Dirac spectrum, the eigenvalues do not 
occur in pairs ±iA, and strictly zero eigenvalues are absent. The Dirac spectrum 
of this theory on the microscopic scale can be analyzed along the same lines as 
the four-dimensional theory. The microscopic spectral density and the Leutwyler- 
Smilga sum rules have been derived [ p.25| , [L26| , the low-energy effective theory 
has been identified p5], pTl], the mass-dependent microscopic spectral density 



has been found [59, 12S], universality has been studied |56|, |129| ], and lattice QCD 



studies have confirmed the theoretical analysis [130]. 



5 UNIVERSAL PROPERTIES OF THE LATTICE QCD DIRAC 
SPECTRUM 

5.1 Unfolding 

In order to compare eigenvalues of any physical system with RMT results, it 



is necessary first to "unfold" the empirical spectrum (see e.g. Ref. [131]). The 
spacing A of the eigenvalues is related to the spectral density by 

A(E) = l/p(E) . (96) 

Unfolding is a local rescaling of the energy scale so that the mean level spacing 
of the unfolded eigenvalues is equal to unity throughout the spectrum. This 
is achieved by splitting the spectral density in an average density, p(E), and a 
fluctuating piece, 

p(E) = p(E) + p a (E) . (97) 

The average spectral density can be obtained in different ways. In some cases, it 
can be computed analytically using semiclassical arguments. However, in most 
cases p(E) is obtained by averaging over many level spacings. There are two 
essentially different procedures to do this, spectral averaging and ensemble aver- 
aging. Spectral averaging is appropriate if only one or a few spectra are available. 
The average level spacing at E is then obtained by averaging over many level spac- 
ings around E. Ensemble averaging is appropriate if there is a large ensemble 
of spectra all drawn from the same statistical distribution. In this case, the av- 
erage spacing can be obtained by averaging the spacing at E over each member 
of the ensemble. In both cases, the separation of the spectral density into an 
average and a fluctuating piece requires a separation of scales that is typically 
achieved only in the thermodynamic limit or in the semiclassical limit. The un- 
folded spectrum, {x n }, is then obtained from the average spectral density (equal 
to the inverse average level spacing) by 

p(X)d\ , (98) 

-oo 

where {E n } is the original sequence of eigenvalues. The difference between the 
two procedures is that spectral averaging yields p(E) separately for each spectrum 
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whereas ensemble averaging yields a single p(E) for all spectra. In both cases, 
the average spectral density of the unfolded spectrum becomes p(x) = 1, so that 
the sequence {x n } has an average level spacing equal to unity. The correlations 
of the levels are always calculated for the unfolded spectrum. 

The two procedures to calculate the average spectral density do not necessarily 
give the same spectral correlations at long distances. The equivalence of the 
two is known as spectral ergodicity. This property has been shown analytically 
for several random matrix ensembles [132, 133]. In general, ensemble averaging 



results in stronger level fluctuations than does spectral averaging [134, 135, 136]. 
In practice a mixture of both methods is often useful. For example, one may 
calculate the average spectral density by ensemble averaging but, in order to get 
better statistics, calculate the correlation functions by both ensemble averaging 
and spectral averaging. Spectral averaging requires that the statistical properties 
of the eigenvalues be stationary over the spectrum, which is not the case for many 
systems and has to be checked each time. 

5.2 Lattice Tests of Chiral Random Matrix Theory 

We have presented a wealth of analytical evidence supporting the statement that 
the local spectral correlations of the Dirac operator are described by universal 
functions. Although the Dirac spectrum is not directly observable in experiments, 
we can compare the predictions of chiral RMT to numerical data obtained by 
lattice gauge simulations. This is the subject of this section. 

This review does not discuss the global spectral density of the QCD Dirac op- 
erator. Most lattice results have large finite-size artifacts. Instanton simulations 
suggest, contrary to random matrix results, that the fermion determinant has a 



significant effect on the global Dirac spectrum [pQ, 122, 137, 135 



5.2.1 BRIEF INTRODUCTION TO LATTICE QCD 

QCD is a renormalizable quantum field theory that must be regularized. This 
can be done by formulating the theory on a discrete lattice with lattice spacing a 



[139]. The largest momentum is then %/a. The quark fields live on the sites and 
the gauge fields on the links of the lattice. If the lattice is finite, one can simulate 
the theory on a computer. This can be done efficiently only in Euclidean space, 
where the gluonic weight function is exp(— Sym) with Sym real. The discretized 
form of the Yang-Mills action Sym is the Wilson action S\y- The full weight 
function of QCD also contains the fermion determinants, which can be expressed 
in terms of the gauge fields. Observables are computed by generating gauge field 
configurations in a Monte Carlo update procedure and averaging an observable 
over many configurations. Since the inclusion of the fermion determinants is very 
time consuming, it is common to use only the gluonic part of the weight function 
in the Monte Carlo updates. This is called the quenched approximation, which 
corresponds to the limit Nf = or, equivalently, to the limit of infinitely heavy 
sea quarks. To make contact with continuum physics, the results of lattice sim- 
ulations must be extrapolated to infinite lattice size (the thermodynamic limit) 
and to zero lattice spacing (the continuum limit). 

Lattice simulations are the main source of nonperturbative information about 
QCD. Unfortunately, the naive discretization of fermions on a lattice leads to the 
so-called doubling problem: The quark propagator has poles at each corner of 
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the Brillouin zone, which gives rise to a total of 2 d species in d dimensions. The 
unwanted 2 d — 1 species can be eliminated by adding to the Dirac operator an 
additional term, the Wilson term, which removes the doublers in the continuum 
limit. However, this term breaks chiral symmetry explicitly. Another possibility 
is the use of staggered (or Kogut-Susskind) fermions where one has only one 
spinor component per lattice site. This maintains a residual chiral symmetry 
but only partially reduces the number of species to 2 d//2 in the a — > limit. 
A no-go theorem by Nielsen and Ninomiya |140| ] states that it is not possible 
simultaneously to solve the doubling problem and have exact chiral symmetry on 
the lattice with a local action. 

Fortunately, there is a way around this theorem: the remnant chiral symmetry 



condition of Ginsparg and Wilson [141] 



^75 + 75^ = 2aV lb RV , (99) 

where R is a spatially local operator that is trivial in Dirac space. In contrast to 
{D, 75} = 0, chiral symmetry is not exact but is recovered only in the continuum 
limit a — > 0. Thus, with Dirac operators satisfying Eq. (|99|), it is possible to get 
rid of the doublers even at finite lattice spacing without violating the Nielsen- 
Ninomiya theorem. Recently, several solutions of Eq. ( |99] ) have been found in the 
overlap formalism 142 1, in the domain- wall formulation in five dimensions [143], 



and in the perfect action approach [144] 



As mentioned in Sec. ggjg , there may be cases in which the anti-unitary sym- 
metries of the various lattice discretizations of the Dirac operator differ from those 
of the continuum operator. In particular, this is true for staggered fermions. In 
SU(2) color, staggered fermions are in the symmetry class of the chSE, whereas 
continuum fermions are in the chOE symmetry class. Staggered fermions in the 
adjoint representation of the gauge group (for any N c ) have the symmetries of the 
chOE, whereas in the continuum limit the symmetries are those of the chSE. The 
Wilson Dirac operator 2?\y does not anticommute with 75 and is therefore not de- 
scribed by any of the chiral ensembles. However, the Hermitian operator 7s£>w 
has the same anti-unitary symmetries as the continuum Dirac operator and is 



described by the corresponding nonchiral ensembles (see also Ref. [145])- Finally, 
Dirac operators obeying the Ginsparg- Wilson condition, Eq. (|99|), have the same 
anti- unitary symmetries as the continuum Dirac operator, but the eigenvalues 



are located on the complex unit circle [62] 



5.2.2 BULK CORRELATIONS 

To measure the bulk spectral correlations, one needs all eigenvalues of the Dirac 
operator, which is represented on the lattice by a finite sparse matrix whose 
dimension is proportional to the lattice volume and to N c . The eigenvalues of 
such a matrix can be obtained using special algorithms, e.g. by the Cullum- 
Willoughby version of the Lanczos algorithm [ |146| ], On the lattice, this was first 
done by Kalkreuter [147|. The numerical effort of this method scales with the 



square of the matrix dimension. For example, in SU(3) on a 10 4 lattice, the 
Dirac matrix has 15,000 distinct positive eigenvalues, which can be computed 
on a typical workstation in about 40 minutes. There are exact sum rules for 
the sum of the squares of all Dirac eigenvalues, which can be used to check the 
numerical accuracy. Because the number of eigenvalues per configuration is very 
large and because the ensemble average can be replaced by a spectral average 
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Figure 3: Number variance S 2 (n), spectral rigidity A^^n), and nearest-neighbor 
spacing distribution P(s) of the lattice Dirac operator with gauge group SU(2). 
Top: Wilson fermions, V = 8 x 12, Nf = 2. Bottom: staggered fermions, 
V = 12 4 , Nf = 4, ma = 0.05. (From Ref. [p9j . Note that in Ref. [p^| , the 
figure of P (s) for Wilson fermions shows data points for staggered fermions. This 
has been corrected in the present figure.) 



under the assumption of spectral ergodicity and stationarity, one needs only a 
few configurations to construct the bulk spectral correlations with great accuracy. 
The bulk correlations have been measured in a number of lattice studies for 

and in instanton liquid 



135, 151, 15 



all three symmetry classes [148, 149, |15 ! 
simulations [137] for (3 = 2. As mentioned above, they are insensitive to the 
number of flavors and the topological charge. In all cases, excellent agreement 
with the predictions of the appropriate random matrix ensemble was obtained 
(see e.g. Fig. ||). The agreement is perfect not only in the strong-coupling regime 
but also at weak coupling. In fact, the bulk spectral correlations are given by 
RMT even in the deconfinement phase [150, 151], indicating that the gauge fields 
retain a sufficient degree of randomness in this region of the phase diagram. 

Spectral ergodicity was investigated in Ref. [135], and the equivalent of a Thou- 
less energy was found for ensemble averaging, whereas spectral averaging resulted 
in complete agreement with RMT correlations over distances as long as several 
hundred average level spacings. 



5.2.3 MICROSCOPIC CORRELATIONS 

Because only the lowest eigenvalues contribute to the microscopic spectral cor- 
relations, a large number of statistically independent spectra are necessary. In 
contrast to the bulk correlations, the microscopic correlations are sensitive to the 
number of flavors and the topological charge. 

The RMT predictions do not contain an energy scale. In order to make com- 
parisons with lattice data, one needs to determine the energy scale 1/V£ to be 
used in the RMT expressions, see e.g. Eq. (p^). This can be done by extracting 
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Figure 4: Valence quark mass dependence of the chiral condensate S(m„) plotted 
as T,(m v )/T, versus m v VY,. The dots and squares represent lattice results by the 
Columbia group ||153| for the values of (3 indicated in the figure. The solid lines 
are chRMT results. (From Ref. fug] .) 



p(0) from the data and applying the Banks-Casher relation, Eq. (pC|). Because 
this procedure makes no reference to RMT, the comparisons of lattice data with 
RMT are parameter-free. 

The first numerical results for the microscopic spectral density p s were obtained 
using instanton liquid configurations for N c = 2,3 and Nf = 0,1,2 [p0[ , and 
the expected agreement with the corresponding random matrix predictions was 
found. The microscopic spectral density was first observed on the lattice via the 
dependence of the chiral condensate on a valence quark mass m v , as studied by 
the Columbia group |153j ]. Figure ^| shows S(m„)/S as a function of m v VT, for 
different values of the coupling constant |116| . The data for different coupling 
strengths in the broken phase fall on a single curve and agree with the chRMT 
result for S(m„) for Nf = \u\ = 0, see Eq. (f)0l). This figure also shows a 
deviation from the chRMT predictions at a scale that is of the order of the 
Thouless energy given in Eq. (|87|). The spectral mass dependence of the chiral 
condensate has recently been studied for different symmetry classes and different 
types of fermions, and a similar quality of agreement has been found [154, |30| . 
The microscopic spectral correlations have also been investigated in great detail 

and the random matrix predic- 



154J, 



on the lattice |5|, [L5|, [TJg, [157], [L5|, [L5£ 
tions were confirmed with very high accuracy for all three symmetry classes. A 
typical example is shown in Fig. ||. 

Although most of these simulations have been performed in the quenched ap- 
proximation, it is also possible to include dynamical fermions. Deviations from 
the Nf = results are observed only if the sea quarks are very light, with masses 
on the scale 1/VT,. Otherwise, the mass term in Eq. (|T|) dominates the small 
eigenvalues in the factors of {i\ n -\- mA. Analytical results for the microscopic 
spectral correlations in the presence of sea quarks with mass ~ 1/VT, were ob- 
tained in Refs. [|160j, 161, 162, 165, 85]. Lattice results in this regime agree very 
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Figure 5: Distribution of the smallest eigenvalue (left) and microscopic spectral 
density (right) of the staggered Dirac operator in quenched SU(2). The dashed 
lines are the predictions of the chSE for Nj = and v = 0. (From Ref. |155(| .) 



well with the corresponding RMT predictions [164]. In the Schwinger model, it 



is numerically feasible to consider massless sea quarks, and the lattice data are 



again well described by RMT [165] 



Recall that the energy scale for the small eigenvalues is 1/VT,. Given the 
agreement of the microscopic spectral quantities with RMT, this means that 



RMT can be used to determine the infinite- volume chiral condensate X [166]. 
This is most easily done by fitting the numerically determined distribution of the 
smallest eigenvalue to the RMT result. Since the lattice volume is known, this 
immediately yields X. 

5.2.4 TOPOLOGY 

The random matrix predictions for the microscopic spectral correlations depend 
on the number of zero modes of the Dirac operator and, thus, on the topological 
charge v of the gauge field configurations. Thus, one should sort the config- 
urations according to their values of v and make the comparison with RMT 
separately in sectors of fixed topological charge. 

The results presented in the previous section were obtained for the staggered 
Dirac operator. If one is interested in topological properties, this operator has 
problems. The zero modes that one would obtain in the continuum limit are 
shifted at finite lattice spacing a by an amount proportional to a 2 [167]. For 



reasonable simulation parameters, this amount is larger than the level spacing 
near zero so that the would-be zero modes are completely mixed with the nonzero 
modes. Thus, one should expect all spectra of the staggered Dirac operator, even 
if computed from gauge field configurations with nonzero topological charge, to 
be described by the RMT results for v = 0. This was indeed observed in the early 
data |116| , |155|| and more recently confirmed in a detailed study (though only in 



strong coupling) [168]. As the continuum limit is approached, the would-be zero 
modes should move toward A = and eventually separate completely from the 
nonzero modes. This effect has recently been observed in the Schwinger model 



in two dimensions [ 169 ] 
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Figure 6: Distribution of the smallest Dirac eigenvalue in the v = 1 sector for 
all three symmetry classes. The data were obtained using the overlap Dirac 
operator on a 4 4 lattice. Solid lines represent the corresponding RMT results. 
(From Ref. gjjg| .) 



Fortunately, Dirac operators satisfying the Ginsparg- Wilson condition (|99| ) do 
not have this problem. The overlap operator [|142f| has exact zero modes even at 
finite lattice spacing, and the microscopic spectral correlations in sectors of fixed 
v are in perfect agreement with the corresponding random matrix predictions 
1 165 , |152| , 154 , |30| , 170 ], Figure [6| shows an example. 



5.2.5 CHIRAL PHASE TRANSITION 



As mentioned in Sec. 3.5, the microscopic correlations are given by chiral RMT 
only if p(0)/V > 0, i.e. if chiral symmetry is spontaneously broken. Chiral sym- 
metry is restored above a critical temperature T c . The way in which p(0)/V 
approaches zero determines the order of the chiral phase transition and, if it is 
second-order, the associated critical exponents. Lattice investigations of these 
questions are very difficult because at finite volume there are no sharp phase 
transitions. Simulations are plagued by finite-size effects, critical slowing down, 
and other problems. It is therefore of great interest to investigate the fate of 
the small Dirac eigenvalues at T = T c analytically in random matrix models, as 



was discussed at the end of Sec. 3.5. There are two recent lattice studies of this 
question [171, 172]. Some of the theoretical expectations were confirmed qualita- 
tively, e.g. the gap in the eigenvalue distribution above T c and the agreement of 
the eigenvalue distribution in this region with the soft-edge predictions of RMT. 
However, further studies are needed to clarify the situation at T = T c . We are 
looking forward to upcoming work in this area. 

5.3 The Thouless Energy and Beyond 

As discussed above, the Dirac spectrum is described by chRMT — or, equiv- 
alently, by the zero-mode approximation of the low-energy effective theory, - 
only below the so-called Thouless energy E c . The theoretical prediction for this 
quantity |a|, pT7|| , 

F 2 
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Figure 7: The data points represent the disconnected susceptibility computed on 
the lattice in quenched SU(3) with staggered fermions (V = 10 4 , (3 = 5.4). The 
solid line is the prediction of chPT, the dashed line the prediction of chRMT. 
(The dashed line is hidden by the data points for u < 10.) (From Ref. |175[| .) 



has been confirmed quantitatively in the instanton liquid model [137, 173] and on 
the lattice [ |174| , [158| , [135| . To compute the Dirac spectrum beyond the Thouless 
energy, the calculation must include the kinetic terms in the low-energy effective 
theory, as discussed in Sec. |4j. The results of such an analysis can again be 
tested by lattice simulations. This has been done for staggered fermions [175], 



Note that in this case, the effective low-energy theory must be modified to take 
into account the symmetries of staggered fermions at finite lattice spacing. A 
convenient quantity to test the theory is the disconnected scalar susceptibility, 
defined by 



disc 



1 / A 1 



X 




M = Jf [ £ {iXk + m){iXl +m y ) ~ Tr ( £ jytt^ ) • ( 101 ) 



This quantity can also be computed in chRMT and in chiral perturbation theory 
(chPT). The result of chRMT is expected to describe the data up to E c . The 
result of chPT should describe the data beyond E c but is expected to break down 
on the scale of the smallest eigenvalue, i.e. for m ~ 1/VTi. These expectations 
are confirmed by the lattice analysis (see Fig. ^ for an example) . 



6 SCHEMATIC MODELS 



We wish to emphasize again that there are two different types of applications of 
RMT to physical problems. RMT may be applied as an exact theory for correla- 
tions of eigenvalues, as discussed in the first part of this review, or it may serve as 
a schematic model for chaos and disorder in physical systems. Below, we intro- 
duce a random matrix model for the chiral phase transition at nonzero chemical 
potential and temperature. Random matrix theories as schematic models are 
mostly used for a description of nonuniversal phenomena. This does not exclude 
that they describe universal fluctuations of the eigenvalues as well. Typically, 
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a random matrix model describes both universal and nonuniversal properties. 
Well-known examples are the Anderson model for localization phenomena [176] 
and the use of random matrix theory in quantum gravity [177, |17~ 



We argued in Sec. 1.1 that chiral symmetry breaking can be understood in 
terms of the stiffness of the Dirac spectrum resulting from interactions of the 
strong color force. Because spectral stiffness is a characteristic feature of RMT, 
it is natural to describe the chiral phase transition in terms of a random matrix 
model. Several different types of schematic random matrix models have been 
introduced. Here we discuss models for the chiral phase transition at nonzero 
temperature T, for the chiral phase transition at nonzero chemical potential (jl, 
and for the phase diagram of QCD. We also discuss random matrix models for 
different types of fermions at u > 0. 

6. 1 Chiral Random Matrix Models for the Chiral Phase Transition at 
Nonzero Temperature 

The original motivation for introducing chiral random matrix models was to ob- 
tain a better understanding of the QCD Dirac spectrum for temperatures around 



the critical temperature for the chiral phase transition [179, 18C[. The idea is 
to split the Dirac operator into the time derivative and a remainder that will be 
replaced by a chRMT, 

V = 7o do + R • (102) 

In a chiral basis with time dependence given by exp[i(2n + 1)itTt], where r 
is the Euclidean time, the first term in T> is diagonal and is given by a direct 
sum of Matsubara frequencies u> n = (2n + \)ttT. The simplest model for T> is 
obtained by replacing the diagonal matrix with positive Matsubara frequencies by 
the identity matrix times an effective Matsubara frequency t, and replacing the 
diagonal matrix with negative Matsubara frequencies by the opposite effective 
Matsubara frequency. After a suitable basis change, this model can be written 
(in the sector of zero topological charge) as |17£, 18C] 



Z Nf (m) = J DWdet N f(V + m)e-^ 2T ^ w ^ , (103) 



where 

iW + it 



v ={ lW i + ti o ) ■ ^ 

We restrict ourselves to the unitary case (j3 = 2) with a complex matrix W 
of dimension N/2. The integral is over the real and imaginary parts of the 
matrix elements of W. For simplicity, we consider only a Gaussian probability 
distribution. The normalization is such that the parameter £ is equal to the 
magnitude of the chiral condensate at zero temperature, 



1 

lim nJ im ^TT dm l °g z N f (m) . (105) 

m^O N-+oo IV ivy t=0 



It can be shown that the effect of the fermion determinant on the macroscopic 
spectral density is subleading in Nj/N, so that p(X) can be calculated in the 
quenched limit. By expanding the resolvent of the Dirac operator, defined by 
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Figure 8: The spectral density of the QCD Dirac operator as a function of an 
effective Matsubara frequency t that models the effect of temperature. (From 
Ref. [I8jj .) 



in a geometric series of the random matrix, it can be shown that g = G/Y> 2 
satisfies a cubic equation |179| , [75], |181| , 

g 3 - 2zg 2 + g{z 2 - t 2 + 1/S 2 ) - z/T? = . (107) 



A variant of the method by which this equation was derived in Refs. |179| , |75j is 
sometimes referred to as the Blue's function method [182]. The spectral density, 
obtained from the discontinuity of the resolvent, is shown in Fig. [8|. We observe 
that chiral symmetry is broken up to t = 1/S with a chiral condensate given by 



G(T, z^0) = 1 - (St) 2 . (108) 

Above this temperature, the spectrum splits into two disconnected regions. Since 
this model has no spacetime dependence, it should not be a surprise that all 
critical exponents are given by their mean field values. From Eq. ( |108[ ), we see 
that P = 1/2. For t 1/S, we find 

G oc z 1 / 3 , (109) 

resulting in another critical exponent 5 = 3. The spectral density at the critical 
point follows by taking the discontinuity of G across the imaginary axis, see 
Eq. ([77]), and is thus given by 

P(A) = — A / . (110) 

The unfolded eigenvalues in this case are given by 

3A \ 4 / 3 m-n 
x k = —X k ' , (ill) 

and a nontrivial scaling limit at t = 1/S is obtained by introducing the micro- 



scopic variable [p8 



z 



AiV 3 / 4 . (112) 
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Starting from the exact analytical result derived in Ref. |7q |, it is possible to take 
the limit N — > oo with the scaling relation (p.l2[) and to obtain an analytical result 



for the microscopic spectral density at the critical point [38, pS[] . We emphasize 
that this result is based on mean field critical exponents and is thus not applicable 
to QCD with nontrivial critical exponents (see also the discussions in Sees. [T^ 
and |5T2^ ). 

As a last application of the partition function ( |103j ) , we mention the explanation 
of the observation |153| ] that the temperature at which chiral symmetry is restored 
is higher for gauge field configurations with a nonzero Z3-phase than for gauge 
field configurations with zero Z3-phase. The explanation in terms of the random 
matrix model ( |103|) is that the lowest Matsubara frequency for the nonzero Z3- 
phases is shifted to a lower value by the phase of the Polyakov loop, and thus the 



critical temperature is higher [181] 



6.2 Chiral Random Matrix Models at Nonzero Chemical Potential 

6.2.1 QCD PARTITION FUNCTION AT NONZERO CHEMICAL POTEN- 
TIAL 

The QCD partition function at nonzero temperature T and chemical potential u 



is given by [183] 



H QCD-M-'V — E a -uN a 

Z(m, u, T) = Tr e r =Y, e T — > ( 113 ) 

a 

where -Hqcd is the Hamiltonian of QCD with eigenvalues E a and N is the quark 
number operator. At zero temperature, only the states with E a /N a < u con- 
tribute to the partition function. We thus expect that the partition function 
is independent of fi below a critical chemical potential \jl c given by the lightest 
baryon mass per unit quark number. Therefore, for [i < u c ~ m^/3 (where m^ 
is the nucleon mass), the baryon density remains zero and the chiral condensate 
is constant. 

The quark chemical potential appears in the Lagrangian in the form fitp'tp = 
^(^70)^ an d is therefore introduced in the Dirac operator by the substitution 

d -> d + (i . (114) 

This substitution destroys the Hermiticity properties of the QCD Dirac operator, 
and the resulting complex phase of the fermion determinant makes Monte Carlo 
simulations impossible. Because of the success of the quenched approximation in 
lattice QCD simulations at u = 0, it is tempting to ignore the fermion determinant 



in this case. However, it was shown [184, |185|] that the critical chemical potential 



for quenched simulations is determined by the pion mass rather than the nucleon 
mass. Obviously, the basic physics of the problem is not visible in the quenched 
approximation. An analytical understanding of this problem was first obtained 
by means of a random matrix model for the Dirac operator of QCD at finite 
density 



6.2.2 A RANDOM MATRIX MODEL FOR QCD AT FINITE DENSITY 

A random matrix model for QCD at nonzero chemical potential model is obtained 
by writing the Dirac operator as ][186|| 



V{n) = a l0 + n (115) 
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and replacing 1Z by a chiral random matrix ensemble. The partition function is 
thus given by 

Z Nf (m) = J DWdet N f(V + m)e-^ 2T ^ w ^ (116) 
with a Dirac operator given by 

iW + n 



D =l.nrt + „ «> >■ (117) 



For QCD with three or more colors (chUE) this Dirac operator has no Hermiticity 
properties. Therefore, the fermion determinant in Eq. ( |116| ) has a complex phase, 
and the partition function cannot be simulated by Monte Carlo methods. 



The partition function with Dirac operator (117) can be rewritten in terms of 
a cr-model. For Nf = 1, the partition function is particularly simple, 

Z(ji) = J dada*[(a + m)(a* + m) - fj, 2 ] n e - n ^ 2 , (118) 

where we have set £ = 1 for ease of notation. For n = N/2 — > oo the integrals 
can be evaluated by a saddle-point approximation. In the chiral limit, we find 



that }186|] 



-_\ for fi>fi c => Z = fi 2n 

a -\ ^l + JP for n < Mc =}► Z = e -"^ 3+1 ) ' ' 

with // c given by the point where the two partition functions are equal, i.e. fx 2 = 
exp(— 1 — ii 2 ), which is solved by /i c ~ 0.53. The vacuum properties of this 
partition function depend on fi for fi < fi c - The reason for this unphysical result 



is discussed below in Sec. 6.2.4. As was demonstrated in Ref. |187[| for (3 = 2, this 



model does not show diquark condensation. However, a random matrix model 



that does show diquark condensation was formulated in Ref. [188] - 

The QCD partition function ( |116|) is a polynomial in m and [i with coefficients 
that can be obtained analytically |189[| . Its properties can be analyzed by means 
of its zeros in the complex chemical potential plane and the complex mass plane. 
Figure ^ shows a first-order phase transition at fj, ~ 0.53. In the chirally restored 
phase, the cut on the imaginary m axis is no longer present. The figure also 
shows the points (stars) where two solutions of the saddle-point equations of the 
cr-model coincide. Indeed, the cuts end on these branch points. 

In lattice QCD at finite density, the zeros of the QCD function in the complex 
\i plane can be studied by means of the Glasgow method [|190|] . This method has 
been analyzed in terms of the above random matrix model J19l[ | with the conclu- 
sion that it requires an exponentially large number of gauge field configurations 
to obtain statistically significant results. 

6.2.3 ZERO TEMPERATURE LIMIT 

The random matrix model discussed in the previous section has attracted a great 
deal of interest, and more sophisticated versions of this model have been proposed 
192] , [189| , [193] , |j"94l , plj . One problem of the model (fTTeJ) is that the chiral 



condensate is /U-dependent below the critical value of the chemical potential. The 
correct zero temperature limit, with a chiral condensate that remains constant 
up to the critical value of the chemical potential, is obtained if all Matsubara 
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Figure 9: Zeros of the partition function in the complex \x and m plane. (From 
Ref. |J|.) 

frequencies are taken into account. We illustrate this with a one-dimensional 
lattice QCD model. 

Writing the Dirac operator as 

Vfa) = 7o(3 + M) + R , (120) 

the /U-dependence of the fermion determinant can be gauged away by means of a 
time-dependent gauge transformation, i.e. 

7o(d + M) + R = ^ [7o9 + R] e» T . (121) 

The /i-dependence is now in the boundary conditions, but for T -> we expect 
that they are not important and that the partition function becomes independent 
of The lesson is that in order to obtain a condensate that is [i- independent, one 
must treat the time derivative exactly, or, in other words, all Matsubara frequen- 
cies have to be taken into account. This was first worked out in detail for lattice 
QCD models at nonzero chemical potential |195(| . Let us discuss the strong cou- 
pling limit of a one-dimensional SU(iV c ) lattice QCD model with Kogut-Susskind 
fermions |196| 1. In that case, the matrix elements of the discretized Dirac operator 
are given by 

V ki = S^Uue" - 5 kH iUie-» , (122) 

where the indices are modulo N, the number of lattice sites. Antiperiodic bound- 
ary conditions result in an extra minus sign for the matrix elements T>in and Djvi . 
The matrices U k i are independent SU(iV c ) matrices. By using gauge invariance, 
one can reduce this partition function to a single SU(iV c ) integral, which can be 
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performed analytically, resulting in the partition function 

sinh[(iV c + l)7Vsinh _1 m 
sinhfiV sinh -1 m 



Z N = 2 cosh(N c N^ + -i . 1 ■ (123) 



The number of lattice sites N can be interpreted as the total number of Matsubara 
frequencies. A sharp transition is obtained in the limit N — > oo. For example, 



the chiral condensate is given by [196] 



v 1 a i rr I N c /y/l + m 2 forsinhu<m, 
km — d m logZ N = i c/ v i (124 

tv^oo N I tor smh fj,> m . 

Similar results have been derived for a large-d strong coupling expansion of lattice 
QCD at finite density |197j 1 (with d the Euclidean dimensionality). The failure 
of the quenched approximation has been reproduced in such lattice models. The 
effect of including all Matsubara frequencies has also been investigated in a ran- 
dom matrix model obtained by replacing the remainder R in Eq. ( |12C| ) by a 
chiral random matrix. After ensemble averaging, one obtains a similar partition 
function with similar conclusions [198]. 



6.2.4 QUENCHING AT NONZERO CHEMICAL POTENTIAL 

The failure of the quenched approximation at nonzero chemical potential was first 
understood analytically in terms of chRMT by Stephanov |186[[ . He showed that 
the quenched limit is the limit Nf — > of a partition function with the absolute 
value of the fermion determinant, 

|det(£>(/x) + m)\ N f , (125) 

rather than 

[det(X>(/i) + m)] N f . (126) 
The absolute value of the fermion determinant can be written as 

det(V + m )det(tf+m*) = det( iW + ^ ,™ ) det ( ~ " .™* ). 

\ m iW'+fiJ \ m iW — fj, ) 

(127) 

Writing the fermion determinant as a Grassmann integral, we observe that the 
quenched partition function can be interpreted as a partition function of quarks 
and conjugate antiquaries. Therefore, in addition to the usual Goldstone modes, 



we have Goldstone modes consisting of a quark and a conjugate antiquark [ 199 , 



186]. Such modes, with the same mass as the usual Goldstone bosons, have a 
nonzero baryon number. The critical chemical potential given by the mass of 
the lightest particle with nonzero baryon number is thus m n /2. The failure of 
the quenched approximation thus has an important benefit: It allows us to write 
down the exact low-energy effective partition function, which is discussed in the 
next section. The product of the two determinants in Eq. ( |127| ) can be written 
as the determinant of a single Hermitian matrix. This procedure is known as 



Hermitization pUO, |20l], pi 



6.2.5 QUENCHED DIRAC SPECTRA 

A chRMT for quenched QCD at nonzero chemical potential is obtained by re- 
placing the determinant in Eq. (116) by its absolute value. This model can be 
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solved analytically in the large- N limit |186| ], However, we will follow a different 
approach by analyzing the corresponding effective theory [202]. The advantage 



of this approach is that the effective theory is as valid for quenched QCD as it is 
for the quenched chiral random matrix model. 

For non-Hermitian matrices, the eigenvalues are scattered in the complex plane. 
Using the fact that 

d z *± = ir5 2 (z), (128) 
z 

where the complex delta function is defined as 5 2 (z) = 5(Re z)5(lm z), we find 
that the two-dimensional spectral density is (up to a normalization constant) 
given by 

p(A) = -d z ,G{z) . (129) 

The resolvent is defined as 

G(z) = J ± v d z logZ, (130) 

and the quenched result is obtained in the limit Nf = of the partition function 
Z = (det 7V /(P(/x) + z)det 7V /(P t (/i) + z*)) . (131) 
The product in Eq. ( |131| ) can be written as 

dflt"/(2?(M) + z)<te»>{tfM + z*) = det ( Wl2N f + Br tWfla £ + ^ ) 

(132) 

with 





B = Bh-B l -\ '■*■< | and C=( ^ z n Nf h ^ 




where we have displayed the degeneracy in flavor space by means of the identity 
matrices I2JV) and 1^,. For z = fj, = 0, the quenched partition function is thus 
invariant under S\Jn(2Nf) x SUi(2iV/). This symmetry is broken spontaneously 
to the diagonal subgroup, SUy(2iVy). At low energies, the effective partition 
function is therefore given by a partition function of Goldstone modes parame- 
terized by matrices U G SU(2iVj). The static effective Lagrangian is obtained 
from the requirement that it should have the same symmetries as the underlying 
microscopic partition function. The mass term was discussed in Sec. [O]. The 
chemical potential term remains invariant under SU(2iVy) x SU(2iVy) transfor- 
mations if at the same time Br and Bl, now considered as independent matrices, 
are transformed as 

Br -> UrBrUk 1 , B L U L B L XJJ} . (134) 

The matrices U in the Goldstone manifold transform as U — > UrIIU^ 1 . To lowest 
nontrivial order in //, we can write two invariant combinations, 

TxUB L U~ l B R and Tr BB . (135) 

The coefficients follow from the conditions that the critical chemical potential is 
equal to one third of the mass of the lightest baryon and that the baryon density 
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should vanish below /i c . For the static limit of the effective partition function, 
we find 

Z= [ DUe -^fTr[U,B][U~\B} + \T,VTr{MU+M^U-^) > ( 136 ) 

Juesv(2N f ) 
with the mass matrix given by 

M=( Z1 P ) • (137) 





\ ' j / 

This partition function can be derived more elegantly by means of a local gauge 



invariance principle fl203| , 204]. For real masses z the Dirac operator satisfies the 
relation det(T)*(fi) + z*) = det(2?(— n) + z), so that the partition function defined 
in Eq. (|131|) and the corresponding low-energy effective partition function given 
in Eq. ( |136|) can be interpreted as the partition function of QCD at finite isospin 
density poll. 

Below the critical chemical potential, u < m n /2, only the vacuum state con- 
tributes to the partition function, so that 

Z = e s ™/(*+**) . (138) 



In terms of the effective theory, the saddle point is at U = 1. The resolvent ( |130|) 
is given by 

G(z) = -l—d z logZ = X. (139) 

JS f V 



Since G(z) does not depend on z*, it follows from Eq. (]12£| ) that the spectral 
density is zero in the region where (|139| ) holds. Because m% = (z + z*)T,/F 2 , see 
Eqs. ( |83|) through (|85|), the condition < m n /2 means that the spectral density 
vanishes everywhere except in a strip 

|Re z\ < J -=— ■ ( 14 0) 

For 2 inside this strip, the Goldstone modes contribute to the partition func- 
tion. In terms of the effective partition function, the saddle point rotates away 
from [7 = 1, leading to a nonvanishing diquark condensate. The rotation of 
the saddle point is a generic feature of low-energy effective partition functions of 



non-Hermitian field theories p0(l 207, 203, 204, 205]. 



These results are in qualitative agreement with spectra of the quenched lattice 



QCD Dirac operator at \x / [184J . The correlations of the complex eigenvalues 



are not so well understood, but the first lattice QCD calculations | 20E ] seem to 



confirm the theoretical expectations [21C, 211, 212, 213 ]. 



6.3 Phase Diagram for the QCD Partition Function 

The random matrix models at nonzero temperature and nonzero chemical po- 
tential introduced in the previous sections can be merged into a single schematic 
chRMT model for the chiral phase transition [|183[] , 



DWdet N f[ ,™ . n iW + iC ) e -^™w< (141) 
■ tW> +iC im I 
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Figure 10: Phase diagram of QCD with two light flavors of mass m as calculated 
from the random matrix model. The almost parallel curves on the wing surface 
are cross sections of this surface with m =constant planes. (From Ref. [|183| .) 



where C is a diagonal matrix with Ckk = t — in. for one half of the diagonal 
elements and C^k = —t — ifi for the other half. In the following, we set £ = 1 for 
simplicity. 

This random matrix partition function can be rewritten as a cr-model. For the 
case Nj = 1 it is given by 



Z(t,/i) 



due 2 



^Trn(o-) 



(142) 



where 



n(a) = aa* - - log[(<r + m)(a* + m)-( f i + it) 2 } 



1 



log[(<7 + m)(a* + m) — (fi — it) 2 ] 



(143) 



For t = this is exactly the a- model discussed in Sec. 3.2.2. At the saddle point, 
the chiral condensate is given by the expectation value of a. For m = the 
saddle-point equation is of fifth order in a, 



a[a 4 - 2(/x 2 



t 2 + ^)a 2 + (v 2 + t 2 ) 2 + v 2 



0. 



(144) 



The critical points occur where one of the solutions of the quartic equation merges 



with the solution a 



0, i.e. along the curve (fi 2 + t 2 ) 2 + n 2 



t l 



tricritical point, three solutions merge. This happens if in addition /i 2 — 1 2 + ^ 



0. At the 
0. 

Figure |10| shows the phase diagram in the /^tm-space. In the m = plane, 
we observe a line of second-order phase transitions and a line of first-order phase 
transitions. They join at the tricritical point. Also joining at the tricritical 
point is a line of second-order phase transitions in the m-direction, which is the 
boundary of the plane of first-order transitions in /itm-space. This line is the 
collection of end points of lines of first-order phase transitions. We expect that 
the critical exponents for such a liquid-gas transition are given by the three- 
dimensional Ising model. The tricritical point was also found in a Nambu-Jona- 



Lasinio model [214, 215]. 
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Because the saddle-point equatio n (|144| ) is of fifth order in a, the critical prop- 
erties of the random matrix model ( |141| ) are very similar to those of a 6 -theory. 
The critical dimension at the tricritical point of such theories is three, so that 
mean field theory, and therefore RMT, describes the correct critical behavior at 
this point. This random matrix model can be considered as the matrix equiv- 
alent of a Landau-Ginzburg functional. The advantage over using the standard 
Landau-Ginzburg theory is that in this case the spectrum of the Dirac operator is 
also accessible. This allows us to study the critical properties of the Dirac eigen- 
values. For example, for /i = 0we have found that the distribution of the smallest 
nonzero eigenvalue of the Dirac operator may serve as an order parameter for the 



chiral transition | 179 |. 



6.4 Random Matrix Triality at \x 7^ 

In previous sections, we have shown that the pattern of chiral symmetry break- 
ing and the correlations of the Dirac eigenvalues are related to the anti-unitary 
symmetries of the Dirac operator. Since the chemical potential occurs only in 
the combination do + /i, the anti- unitary symmetries at /j, 7= are the same as 
for zero chemical potential. Thus for QCD with two colors and fermions in the 
fundamental representation the Dirac operator is real (/3 = 1), whereas for QCD 
with adjoint fermions the Dirac operator is quaternion real {(5 = 4). In the first 



case, the Dirac operator has the structure [ 195 ] 



V=[ _° „ W y ) with W real. (llo) 



-W 1 + n 
In the second case, T> is given by 

V = ( _|_ ^ ^ Q ^ J with W quaternion real. (146) 

In this case, the quaternion real matrix elements of W satisfy the reality relation 
Wj^i = cx^Wki^vt- I n both cases, it is easy to show that the fermion determinant 
is real. Furthermore, because of 

&et N fW&et N fW T = &et 2N fW for (3 = 1, (147) 
det^ W det N f W ] = &et 2N f W for (3 = 4, (148) 

the flavor symmetry group is enhanced to \J{2Nf) in both cases (for v = O)]] 
The full flavor symmetry is broken spontaneously to Sp(2iVy) and 0(2Nf), re- 
spectively, and is broken in the same way by the mass term. The chemical 
potential breaks the symmetry according to U(2iV/) — ► U(iVy) x U(iVj) in both 
cases. For both (3 = 1 and /3 = 4 the pseudoreality of the Dirac operator leads 
to Goldstone modes with a nonzero baryon number in the same way we have 
seen for the phase quenched partition function. An exact low-energy effective 
partition function valid to lowest order in m and fi 2 can be written down [204]. 



3 For p = 4 and fi = 0, the Dirac operator satisfies the relation (CV) T = —CD whereas 
(C^7o) T = Cfijo (with C the charge conjugation matrix). Therefore, for j3 = 4 and /i / 0, 
the Pfaffian of the Dirac operator is not defined and the fermion determinant Aet N{ ^ 2 {V + fi-yo) 
cannot be expres sed a s an integral over Nf Majorana fermions with flavor symmetry U(iV/), as 
discussed in Sec. |2.2.S| , Instead, we write det Nf (T> + A*7o) as an integral over Nf Dirac fermions 
with flavor symmetry U(2iV/). 
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Figure 11: Scatter plot of the real (x) and imaginary (y) parts of the eigenvalues 
of the random matrix Dirac operator at nonzero chemical potential. The val- 
ues of (3 and \x are given in the figures. Full curves show the analytical result 
for the boundary. The line along the imaginary axis for (3 = 1 represents the 



accumulation of eigenvalues on this axis. (From Ref. |193| .) 



For an elaborate discussion of the symmetries, their breaking pattern, and the 
low-energy effective theory in both cases, see Ref. [204]. 

For an even number of flavors, dynamical Monte Carlo simulations are possible 
for (3 = 1 and (3 = 4, though not for [3 = 2. However, quenched simulations 
have been performed for all three classes. A cut along the imaginary axis below 
a cloud of eigenvalues was found in instanton liquid simulations [ 21C ] for N c = 2 
at u. ^= 0, which corresponds to (3 = 1. In lattice QCD simulations with staggered 
fermions for N c = 2 [217], a depletion of eigenvalues along the imaginary axis 
was observed, whereas for N c = 3 the eigenvalue distribution did not show any 
pronounced features [184]. 

A chiral random matrix model for the Dirac operator at [i ^ in each sym- 
metry class is obtained by drawing the matrix elements of W from a Gaussian 
probability distribution. In the quenched approximation, the spectral proper- 
ties of the random matrix Dirac operator of Eq. ( |117| ) can easily be studied 
numerically by diagonalizing a set of matrices with the probability distribution 
of Eq. (|116| ). Figure 11 shows results [193] for the eigenvalues of a few 100 x 100 
matrices for jj, = 0.15 (dots). The solid curves represent the analytical result 
for the boundary of the domain of eigenvalues derived in Ref. |186(| for (3 = 2. 
However, the analysis can be extended [193] to (3 = 1 and (3 = 4, and with the 
proper scale factors, the solution is identical. 

For (3 = 1 and (3 = 4 we observe exactly the same structure as in the previously 
mentioned (quenched) QCD simulations. There is an accumulation of eigenvalues 
on the imaginary axis for (3=1 and a depletion of eigenvalues along this axis for 
(3 = 4:. The depletion can be understood as follows. For fj, = all eigenvalues 
are doubly degenerate. This degeneracy is broken at /i ^ 0, which produces the 
observed repulsion between the eigenvalues. 

The number of purely imaginary eigenvalues for (3 = 1 scales as y/~N and is 
thus not visible in a leading-order saddle-point analysis. Such a y/N scaling is 
typical for the regime of weak non-Hermiticity first identified by Fyodorov et 
al. [211]. Using the super symmetric method of RMT, Efetov f]2Q7f obtained 
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the vN dependence analytically. The case = 4 was also analyzed analytically 
[p!8| , with results that are in complete agreement with the numerical simulations. 
Obviously, more work is needed in order to arrive at a complete characterization 



of universal features [212] of the spectrum of non-Hermitian operators. 



7 RELATED MODELS AND RELATIONS TO OTHER FIELDS 

Random matrix theory has been used extensively in many different fields includ- 
ing nuclear physics, atomic and molecular physics, condensed matter physics, 
quantum chaos, quantum gravity, and mathematical physics. Before reviewing a 
few applications that might have some relation to QCD, we briefly discuss another 
class of matrix models for the QCD partition function. For a more comprehensive 
review of the material in this section, see Refs. 



18, 21S, 22C 



7.1 One-Plaquette Models of Lattice QCD 

One class of matrix models for QCD are lattice QCD partition functions on a 2 d 
lattice. The simplest model in this class is the one-plaquette partition function 
for pure gauge theory in d = 2. This model is known as the Brezin-Gross-Witten 
model p23[| and is defined by 



Z(J, Jt) 



DUe 



Tr(JC/t + jtf7) 



t/elX(7V) 



(149) 



where the integral is over the Haar measure of U(iV) and J is an arbitrary complex 
source term. We have shown that this partition function can be expressed as a 
determinant of modified Bessel functions |m], 23, j93"fl . For J a multiple of the 
identity (with the interpretation of J ~ 1 /g 2 , where g is the Yang-Mills coupling 
constant), this model was solved analytically in the large- iV limit by Gross and 
Witten. They found a third-order phase transition at a critical value of J. A 
large- N chiral phase transition is obtained by extending this model with a fermion 



determinant [224]. Despite considerable effort, the generalization of this model 



to more than two dimensions, which is known as the Eguchi-Kawai model [225], 
could not be solved analytically [p26|. 

The Brezin-Gross-Witten model [223] can be rewritten as a generalized Kont- 
sevich model, which has received a great deal of attention in the theory of exactly 
integrable systems [227, 222]. There are many other relations between RMT and 
the theory of exactly solvable systems. Among others, the asymptotic properties 



of correlation functions can be derived by means of conformal field theory [228] 
An interesting overview in the context of string theory is given in Ref. [227]. 



Another class of related models is the Kazakov-Migdal model [22£], which is 
also known as induced QCD. In its simplest form it is defined by 



(150) 



where the H (x) are Hermitian matrices and the U(x) are unitary matrices. The 
sum is over a lattice in d dimensions. A zero-dimensional form of this model was 
proposed as a model for the transition between Poisson statistics and Wigner- 
Dyson statistics [230, 173]. 
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7.2 Universal Conductance Fluctuations in Disordered Mesoscopic 
Systems 

In condensed matter physics, a mesoscopic system is a system whose linear size is 
larger than the elastic mean free path of the electrons but smaller than the phase 
coherence length, which is essentially the inelastic mean free path. A typical 
size is about 1 /im. The conductance g of mesoscopic samples is closely related 
to their spectral properties. Using a scaling block picture, Thouless found that 
in the diffusive regime, the conductance is given by g = Eq/A, where Ec/h is 
the inverse diffusion time of an electron through the sample and A is the mean 



level spacing [231]. This can be rewritten as g = (N(Ec)), where {N(E)) is the 
mean level number in an energy interval E. Thus the variance, (Sg 2 ), of the 
conductance is related to the number variance, S 2 , of the energy levels |232| 1. 

Low-temperature experiments have been performed in which the conductance 
of mesoscopic wires was measured as a function of an external magnetic field. The 
observed fluctuations in g are of the order of e 2 /h, independent of the details of 
the system (shape, material, etc). These are the so-called universal conductance 
fluctuations ]233f . One can understand this phenomenon qualitatively by esti- 
mating the number fluctuations of the electron levels using RMT results. Both 
the magnitude of the fluctuations and their universality can be obtained through 
the transfer matrix method. An interesting numerical result is that the density 
of eigenvalues of the transmission matrix in the Hofstadter model for universal 
conductance fluctuations can be described in terms of the microscopic spectral 



density of the chUE [234, E3 



7.3 Anderson Localization 

In more than two dimensions, a good conductor becomes an insulator when the 
disorder becomes sufficiently strong. This phenomenon is called Anderson local- 
ization. In the localized phase, the wave function of the electron is not described 
by Bloch waves, but by a localized form that decays exponentially, 

iff(r) ~ e~ r/Lc . (151) 

The length scale L c is known as the localization length. This phenomenon was 
first described by the Anderson model |176| ] , which is a hopping model with a 
random potential on each lattice point. The dimensionality of the lattice plays an 
important role. It has been shown that in one dimension all states are localized. 
The critical dimension is two, whereas in three dimensions there is a derealization 
transition at an energy El. The states below El are localized whereas the states 
above El are extended, i.e. with a wave function that scales with the size of 
the system. The eigenvalues of the localized states are not correlated, and their 
correlations are described by the Poisson distribution. 

An interesting question is whether Dirac eigenfunctions can be exponentially 
localized. Parisi argued that localized states can only occur in quenched systems 



1 236]. In QCD this argument goes as follows. If the eigenfunctions were spatially 
localized, the eigenvalues would be uncorrelated, and there would be no repulsion 
between the eigenvalues. Because of the fermion determinant and the measure 
in Eq. (|34]), the eigenvalues would be repelled from the origin. Since there would 
be no mechanism to compensate for this repulsion, p(0)/V and the chiral con- 
densate would be zero. Therefore, if chiral symmetry is spontaneously broken, 
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the eigenfunctions of the Dirac operator must be spatially extended. Indeed, this 
has been found for the wave functions of the Dirac operator for a gauge field 
ensemble given by a liquid of instantons |137| ] as well as in lattice QCD [ |237| . 

In the extended domain, the situation is more complicated. An important 
energy scale is the Thouless energy ||231|| , which is related to the diffusion time 
of an electron through the sample (see Sec. 7.2). With that diffusion time given 



by L /D (the diffusion constant is denoted by D), the Thouless energy is [23S 



E c = ^ ■ (152) 

On time scales larger than h/E c , an initially localized wave packet diffuses all over 
phase space. If this wave function ip(t) att > h/E c is expressed as a superposition 
of eigenfunctions 4>i(t = 0), many of the overlaps {ip(t)\4>i) are nonzero. Therefore 
we expect that for energy differences below E c the eigenvalues are correlated 
according to RMT. 

The diffusive behavior of electrons is described by Goldstone modes, called 



diffusons. Classically, they satisfy a diffusion equation |238| , 239 , 102]. The 
diffusion modes can be described in terms of a cr-model, which has the same 
structure as the low-energy effective theory for the QCD partition function. The 
Thouless energy corresponds to the energy scale below which the fluctuations of 
the zero-momentum modes dominate the fluctuations of the nonzero-momentum 



modes, i.e. the scale discussed in Sec. 4.4, below which the QCD Dirac spectra 



are given by chRMT. The analysis of classically chaotic quantum systems has 
made it clear that spectra are described by RMT if and only if the corresponding 



classical motion is chaotic | ]240| , 241]. If we interpret the Dirac operator as a 
Hamiltonian in four Euclidean dimensions plus one artificial time dimension, 
we thus conclude that the classical motion of the quarks is chaotic^ [^] and 
therefore diffusive in four Euclidean dimensions and one artificial time dimension, 
in agreement with the interpretation of Goldstone modes as diffusion modes [{2391, 



243]. The equivalent of the diffusion constant can be read off immediately by 



comparing the two a- models, or equivalently by comparing the scale of Eq. (p7\ ) 
with the Thouless energy, and is given by D ~ F 2 /£ [137, [244|| . Quantum 



mechanically, the Thouless energy can be interpreted as the spreading width of 



the exact eigenstates over the states of the noninteracting Hamiltonian 245 , 246]. 



7.4 Non-Hermitian Models with Disorder 

In the past few years, several non-Hermitian models with disorder have been 
introduced in the literature. We mention only the simplest model, which was 



motivated by the study of flux-line pinning in superconductors |247], |24q | . Its 
Hamiltonian is given by 

H = ^-(p + ih) 2 + V(r), (153) 
2m 

where V(r) is a random disorder potential, p is the momentum operator, and ih 
is a constant imaginary vector potential. The new feature of such models is that 
a localization-delocalization transition can occur even in one and two dimensions. 
Although the wave functions corresponding to real eigenvalues remain localized, 

4 There is an elabor ate literature on the classical dynamics of gauge fields. For a discussion 
of this topic, see Ref. |242f . 
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the wave functions of complex eigenvalues can be extended. For a detailed dis- 
cussion of such models, see Refs. [201, 249 1 . 



The dependence on h in the Hamiltonian H occurs only in the combination 
p + ih, i.e. in exactly the same way as the chemical potential in the QCD Dirac 
operator. There are two important differences from QCD: (a) the operator H 
does not have a chiral structure, and (b) the disorder is uncorrelated Gaussian and 
quenched. The connection between derealization in the Hatano-Nelson model 
( |153j ) and diquark condensation in quenched QCD at nonzero chemical potential 
is not yet understood and deserves further attention. 



After the initial work of Ginibre [210], non-Hermitian RMT has received a 



great deal of attention in the mathematics literature as well. Probably the best 



overview of results in this area is in the book by Girko [ 200 1 . 



7.5 Andreev Scattering 

The term Andreev scattering (or Andreev reflection) refers to a process in which 
an electron hits the interface between a normal metal and a superconductor and 



is reflected as a hole with the opposite momentum (or vice versa) | 250 |. Stated 
differently, two electrons can tunnel through the interface so that a Cooper pair 
is added to the superconducting condensate (or removed from it). This process 
can be described in a microscopic mean field model by the Bogoliubov-deGennes 
Hamiltonian, which can be written in matrix form as 

_ B a t), (154) 

where A (— A T ) represents the Hamiltonian for particles (holes) and B represents 
the pairing field. The requirement that H be Hermitian means that A must be 
Hermitian. If the system under consideration is invariant under time reversal 
or spin rotation or both, there may be additional symmetries. If the system is 
invariant under spin rotation, B is symmetric (in this case we consider the Hamil- 
tonian in the space of spin-up states only). If, in addition, the system is invariant 
under time reversal, the elements of H are real. Otherwise, they are complex. On 
the other hand, if the system is not invariant under spin rotation, B is antisym- 
metric. In this case, if the system is invariant under time reversal, the elements 
of H are quaternion real. Otherwise, they are complex. This classification es- 
tablishes four new random matrix ensembles p51| , |33|| . The microscopic spectral 
correlations of these ensembles are identical to those of the chiral ensembles if 
the parameters (the Dyson index (3 and the number of massless flavors Nf) are 
chosen appropriately. 

7.6 Mathematical Physics and Quantum Gravity 

RMT has received a great deal of attention as a problem in mathematical physics. 
We mention the relation between universal behavior in RMT and the asymptotic 



behavior of orthogonal polynomials |42j, |43|, 55, |252j, [253||, the relation between 



the classification of random matrix ensembles and the Cartan classification of 



symmetric spaces [34], and the theory of Riemannian superspaces j34[. Vari- 
ous methods for the solution of random matrix problems have been proposed. 
We mention the mapping to a gas of noninteracting fermions, the Coulomb gas 



method [27|, the Brownian motion method [27], the replica method [107], the 
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orthogonal polynomial method [17|, the super symmetric method |102|, and the 
operator method |177| ]. 

Generally, the problem for (3 = 2 is much simpler than for = 1 and (3 = 4. 
Nevertheless, a great deal of progress has been made for (3 = 1 and (3 = 4. 
We mention relations between the kernels of the correlation functions |47], 48, 
|S2] , 83] and the relation between massless and massive correlators [85]. Novel 
mathematical methods have been developed for (3 = 1 and (3 = 4. We mention 

and the extension 



255, 43, 



only the skew-orthogonal polynomial method [25 
of an operator method for (3 



2 fl81|, £56|, 

There exists an elaborate literature on the random matrix formulation of 2- 
d quantum gravity (see e.g. the recent reviews by Abdalla et al. [175] and Di 
Franceso et al. |177f| ). The idea is that the partition function, which is a sum 
over all metrics, can be rewritten as a sum over random surfaces. The sum over 
the dual graphs of a discretization of these random surfaces can be rewritten in 
terms of a random matrix partition function that can be analyzed with standard 
random matrix methods. 



8 CONCLUSIONS 

Chiral random matrix theory is the archetypal model for the spontaneous break- 
ing of chiral symmetry. In chRMT, chiral symmetry is spontaneously broken for 
any finite number of massless flavors, whereas in QCD, breaking is believed to 
occur only below a certain number of massless flavors, perhaps as few as four. In 
this review, we have studied chiral symmetry breaking from the perspective of the 
eigenvalues of the QCD Dirac operator, where broken chiral symmetry implies 
that the smallest eigenvalues are spaced as 1/V . The robustness of chiral sym- 
metry breaking in chRMT can then be understood as the crystallization of the 
eigenvalues by the long-range random interactions. The number of eigenvalues 
in a sequence containing TV eigenvalues fluctuates on average only by \/IogN /ir 
as opposed to y/~N for uncorrelated eigenvalues. For uncorrelated eigenvalues, we 
do not expect chiral symmetry breaking for any number of massless flavors. The 
conclusion is that chiral symmetry breaking in QCD requires strongly correlated 
Dirac eigenvalues but not quite as strongly correlated as in chRMT. To make 
this conclusion more quantitative, we have formulated an effective theory for the 
QCD Dirac spectrum, which, in addition to the usual Goldstone modes, contains 
Goldstone bosons of spectral quarks that are "ghost" quarks introduced to probe 
the Dirac spectrum. Like the usual chiral Lagrangian, this Lagrangian consists 
of a kinetic term and a mass term. An important scale is the spectral mass for 
which the long-wavelength fluctuations of these two terms are of equal order of 
magnitude. In the theory of disordered condensed matter systems, this scale is 
known as the Thouless energy. For spectral masses below the Thouless energy, 
the mass dependence is given by the zero-momentum sector of the effective the- 
ory, and it coincides with that of chRMT in the limit of large matrices. The 
deviations from chRMT, given by the contributions of the nonzero-momentum 
modes, increase the fluctuations of the eigenvalues. For an increasing number 
of flavors, the position of the smallest eigenvalues moves to the right, whereas 
the slope of the spectral density increases with Nf. However, chiral symmetry 
remains broken for any value of Nf . The chiral condensate is simply a parameter 
of the effective partition function. 
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The statistical properties of the QCD Dirac eigenvalues have been investigated 
by numerous lattice QCD simulations, and the chRMT predictions have been 
verified in great detail. The behavior of Dirac spectra in the domain beyond the 
Thouless energy has been less well studied, but in all cases, the predictions of the 
effective theory for the Dirac spectrum agree well with lattice simulations. 

Inspired by the successes of chRMT as an exact description of the statistical 
properties of the QCD Dirac spectrum, we have constructed a schematic model 
for the chiral phase transition. Although this approach does not provide rigorous 
results, it has been useful in advancing a qualitative understanding of the chiral 
phase transition at both zero and nonzero chemical potential. Examples are the 
properties of the Dirac spectrum in the approach to the critical temperature, the 
failure of the quenched approximation at fj, > 0, and the phase diagram of the 
QCD partition function in the \xTm plane. 

Finally, RMT encapsulates a duality between order and chaos. Complete mix- 
ing of all degrees of freedom can be described by a single effective degree of 
freedom. This is the progress we have made toward an understanding of the 
complexity of the QCD vacuum. 
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